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ABSTRACT 

We present an analysis of the s-process elements Zr, Ba, and La and the r-process element Eu in a 
sample of 50 stars in 19 open clusters. Stellar abundances of each element are based on measures of a 
minimum of two lines per species via both equivalent width and spectrum synthesis techniques. We 
investigate cluster mean neutron capture abundance trends as a function of cluster age and location 
in the Milky Way disk and compare them to results found in other studies in the literature. We find 
a statistically significant trend of increasing cluster [Ba/Fe] as a function of decreasing cluster age, in 
agreement with recent findings for other open cluster samples, supporting the increased importance of 
low mass AGB stars to the generation of s-process elements. However, the other s-process elements, 
[La/Fe] and [Zr/Fe] do not show similar dependences, in contrast to theoretical expectations and the 
limited observational data from other studies. Conversely cluster [Eu/Fe] ratios show a slight increase 
with increasing cluster age, although with marginal statistical significance. Ratios of [s/r]-process 
abundances, [Ba/Eu] and [La/Eu], however, show more clearly the increasing efficiency of s-process 
relative to r-process enrichment in open cluster chemical evolution, with significant increases among 
younger clusters. Lastly, cluster neutron capture element abundances appear to be independent of 
Galactocentric distance. We conclude that a homogeneous analysis of a larger sample of open clusters 
is needed to resolve the apparent discrepant conclusions between different studies regarding s-process 
element abundance trends with age to better inform models of galactic chemical evolution. 
Subject headings: Galaxy: abundances — stars: abundances 



1. INTRODUCTION 

It is well established that open clust ers do not 
exhibit an age-metallici ty relation (e.g., iFriell I1995L 
iCarrera &: Panciiiol 120 111 ) . The fact that old, super- 
solar metallicity clusters (e.g., NGC 188, NGC 6791) and 
young, sub-solar metallicity clusters (e.g., NGC 2324) 
can be found in the disk seems to indicate that the loca- 
tion of a cluster's birth is more important than when it 
formed in terms of the signat ures of disk chemic al evo- 
lution imprinted on its stars (Salaris et al.|[2004T ). This 
is seen in the radial metallicity gradient in the disk, as 
trace d by open clusters , cepheids, and other populations 
(e.g., lYong et al.ll2012f ). The distribution of abundances 
of other chemical elements (such as the alpha-, light and 
Fe-peak elements) also appear to trace that of iron very 
well, with open clusters having [X/Fe] ratios that are es- 
sentially indepen dent of cluster age and galactocentr ic 
radius (R gc ; e.g.. iMagrini et all [20091 iFriel et al.ll2010T Fl. 
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The findings of studies of open cluster neutron-capture 
element abundances, however, stand in exciting contrast 
to the behavior seen for the other elements. Neutron 
capture species, generated in low mass AGB stars and 
massive stars via the strong and weak s-processes and 
(likely) in type II supernovae via the r-process (e.g. , 
I Busso et al.lfl99l IBusso et all [2001 iSneden et aDl2008l 
iPignatari et al.ll2010f ). have only recently been studied 
in decently large, homogeneous samples of open clusters, 
and are further useful probes of the mass range of stars 
that contributed to the chemical evolution of the disk. 
iD'Orazi et al.l (|2009t ) reported an unexpected trend of in- 
creasing cluster [Ba/Fe] with decreasing cluster age that 
cannot be explained by chemical evolution models using 
standard AGB star s-process yields. An analysis of the 
s-process elements Y , Ce, Zr and La in open clusters by 
iMaiorca et al.l (|2011| ) found an analogous trend for these 
elements, though it should be noted that the result is 
most robust for Y and Ce. Given that this trend is seen 
for both light s-elements (Y, Zr) and heavy s-process el- 
ements (Ba, La, Ce) in these studies, one can conclude 
that the cause of this increased enrichment in younger 
stars affects all elements produced by the main s-process 
in the same way. 

in galaxy disks via various mechanisms c an greatly in fluence radial 
metallicity gradients (e.g., Roskar et al. (2008) and Locbman et al. 
(2011)), to our knowledge it has not been demonstrated that open 
clusters are subject to radial migration. It goes without saying the 
usefulness of cluster element abundance trends with R gc as probes 
of disk chemical evolution hinges on the answer to this question. 
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Table 1 

Cluster Sample 





Age a 


d 


R b 


Ref. c 


No. of Stars 


Telescope 


Cluster 


(Gyr) 


(kpc) 


(kpc) 










Be 17 


10.1 


2.7 


11.7 


F05 


3 


KPNO 


4m 


Be 18 


5.7 


5.4 


13.7 


Y12 


2 


Keck 




Be 21 


2.2 


6.2 


14.7 


Y12 


2 


Keck 




Be 22 


1.3 


6.2 


14.4 


Y12 


2 


Keck 




Be 32 


5.9 


3.2 


11.6 


FJP10 


2, 2 


KPNO 


4m, Keck 


Be 39 


7.0 


1.3 


11.9 


FJP10 


1 


KPNO 


4m 


M 67 


4.30 


0.8 


9.1 


FJP10 


3 


KPNO 


4m 


JNGO 188 


6.3 


1.7 


9.4 


FJP10 


1 


KPNO 


1m 


NGC 1193 


4.2 


5.8 


13.6 


FJP10 


1 


HET 




NGC 1245 


1.1 


3.0 


11.1 


Jll 


1 


KPNO 


4m 


NGC 1817 


1.1 


1.5 


10.0 


J09 


2 


KPNO 


4m 


NGC 1883 d 


0.7 


3.9 


12.3 


J09 


2 


KPNO 


4m 


NGC 2141 


2.1 


3.9 


12.2 


J09 


2 


KPNO 


4m 


NGC 2158 


1.9 


4.0 


12.8 


J09 


1 


KPNO 


4m 


NGC 2194 


0.9 


1.9 


10.3 


Jll 


2 


KPNO 


4m 


NGC 2355 


0.8 


1.9 


10.3 


Jll 


3 


KPNO 


4m 


NGC 6939 


1.2 


1.8 


8.4 


A04 


1 


KPNO 


4m 


NGC 7142 


4.0 


1.9 


9.2 


J08 


1 


KPNO 


4m 


PWM 4 


7.0 


7.2 


13.3 


Y12 


1 


Keck 





a Adopted from lSalaris et al.1 £2004) 
b R R = 8.5 kpc 



c Refer ences for dis t ances : A04: [Andrcuzzi ct al. (2004]}; F05 = IFriel et aHl2005T); 
FJP = IFriel et al.l 120101); J08 = IJacobson et al.l ffixM ) ; J09 = IJacobson et al.l 
l|2009l 1; Jll = IJacobson et al.l 1)20111 ): Y12 = I Yong et*aU f20T2T ) 
d N GC 1883 is not include d in the sample of Salaris ct aQ. This value is taken 
from lVillanova et al.1 ( 120071 ) (uncertainty 0.07 Gyr). 



D'Orazi et al. interpreted this anticorrclation with age 
as an indication that low-mass (<1.5 M Q ) AGB stars 
contributed more to the chemical evolution of the disk 
than theoretical models previousl y indicated. This sce - 
nario has been further explored by lMaiorca "eTaII(pr)T2T ). 
who presented revised s-process production models for 
low-mass stars in which the neutron source was larger 
by a factor of four than in previous models. Applying 
the s-process yields produced by this model in a chemi- 
cal evolution study, they could reasonably reproduce the 
observed trend of [s/Fe] with age seen by open clusters. 

It is very intriguing that a recent study of solar-type 
stars in the solar neighborhood (for which there are Hip- 
parcos distances and therefore, isochrone ages) has also 
found trends of increasing [s/Fe] with d ecreasing age. 
Indeed , the dwarf stars in the study by Ida Silva et al.l 
(2012) show a very similar pattern of [s/Fe] with age to 
the open clusters in the D'Orazi et al. and Maiorca et al. 
studies, most notably for the elements Ba and Zr, and 
especially for stars younger than the age of the Sun (see 
their Figure 11). That said, the resul ts of another open 
cluster study show conflicting results. lYong et all (|2012[ ) 
presented a homogeneous chemical abundance study of 
11 open clusters, and in combination of results from the 
literature, investigated trends of cluster element abun- 
dances as a function of R gc and cluster ag^EL While they 
found trends of increasing [Ba/Fe] and [Zr/Fe] with de- 
creasing cluster age, they were not as strong as those 
found by D'Orazi et al. and Maiorca et al. (hereafter 

5 It should be noted for the present discussion that only a sub- 
set of the literature sample they examined had neutron capture 
element abundances reported, a large fraction of which were from 
their own work (see their Table 13), and so the usual risks inher- 
ent in interpreting trends in inhomogeneous datasets are not so 
important here. 



D09 and iMli!)- What is more, they found a trend of in- 
creasing [La/Fe] with incre asing cluster age, in the direct 
opposite sense of the lMllI result. 

Given the important role open clusters perform in 
tracing the chemical evolution of the disk, it is imper- 
ative to firmly establish the existence and behavior of 
their element abundance trends with age. To that end, 
we present here an analysis of neutron capture element 
abundances in a sample of 19 open clusters spanning in 
age from ^700 Myr to 10 Gyr. Our study has many stars 
in co mmon with those of Yong et al. (2012; hereafter 
lY12h , but we have undertaken a completely independent 
analysis of them. Furthermore, only two clust ers in our 
sample are common to the studies of lD09l and lMlll . and 
therefore this study is an independent investigation of 
cluster abundance trends with age. The elements inves- 
tigated in this study are primarily s-process elements: Zr, 
Ba and La. We also report abundances for the r-process 
element europium (Eu). 

2. SAMPLE 

The cluster sample used in this work is shown in Ta- 
ble Q] Much of our stellar sample has been presented in 
previous work, and is comprised of echelle s pectra from 
KPNO 4m and Hobby-Eber l y tele s copes (see IFriel et al.l 
[20051 [2010L IJacobson etall [20081 I2009T ). This sample 
includes confirmed radial velocity cluster member red 
clump and giant stars in the clusters Be 17 (3 stars), 
Be 32 (2 stars), Be 39 (4 stars), M 67 (3 stars), NGC 
188 (4 stars), NGC 1193 (1 star), NGC 1817 (2 stars), 
NGC 1883 (2 stars), NGC 2141 (2 stan0), NGC 2158 (1 

6 The spectrum of star 1348 was obtained by Yong et al. 2005 
using a similar spectrograph setup to the data described here. This 
spectrum was kindly given to us by D. Yong (2008, pr ivate com- 
munication) and was analyzed in Jacobson ct al. (2009). 
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Figure 1. Sample portions of the spectra for stars 1007 (dashed line) and 1348 (solid line) of NGC 2141. The two stars have identical 
atmospheric parameters and metallicities, but 1348 shows enhanced Zr and La abundances relative to 1007, although their Ba and Eu line 
strengths are very similar. 



star) and NGC 7142 (4 stars). The reader is referred to 
references above for details regarding the observations, 
reductions, and analysis of these data and information 
about the stellar targets. 

Additional high resolution spectra of stars in clusters 
Be 18, Be 21, Be 22, Be 32, and PWM 4 obtained with 
the HIRES spectrograph on Keck were kindly given to us 
by D. Yong (two stars in each cluster save one in PWM 
4). These outer disk clusters serve as an important sup- 
plement to our sample, and a co mparison of the analysis 
here to that presented by IY12I provides important in- 
formation about systematic differences between different 
cluster samples in the literature along the lines investi- 
gated by us previously (see Appendix in lFriel et al.H2010l 
hereafter FJP10). We refer the reader to lY12l for details 
regarding the observations and reductions of these spec- 
tra. 

We obtained new observations of stars in open clusters 
NGC 1245, NGC 2194, NGC 2355 and NGC 6939 using 
the echelle spectrograph on the KPNO 4m Mayall tele- 
scope in 2010 November. Cluster radial velocity mem- 
bers were selected for high resolution m ulti-order spec- 
trosco pic follow-up from the sample of Jacobso n et al.l 
(|2011| ). These data are similar to those obtained by us 
in previous runs (e.g., FJP10): the long red camera along 
with the 316 line mm -1 grating and a 1" slit resulted in 
spectra with resolving power R(A/AA) ~ 28,000. The 



full wavelength range is 5000-8100 A across 24 orders. 
Typically 3-6 exposures of no more than 2700 s were ob- 
tained for each target to minimize the impact of cosmic 
rays. In addition, high signal-to- noise (S/N) spectra of 
radial velocity standard stars and early type stars (used 
for subtraction of telluric lines) were obtained each night, 
along with the standard bias, dome fiat, and arc lamp 
calibration frames. 

These data were reduced using standard IRAI0 pro- 
cedures, as described in FJP10. After bias subtraction, 
scattered light removal and flat field division, individ- 
ual data frames were process ed through "L. A Cos- 
mic"0 to remove cosmic rays (jvan Dokkuml l200ll spec- 
troscopic version) . Individual orders were then extracted 
and converted to one-dimensional spectra and disper- 
sion corrected. Individual object spectra were combined 
to form a single, high S/N spectrum for each star and 
then continuum- normalized. Signal-to-noise (S/N) ra- 
tios of the summed spectra range from 60 to 120 per 
pixel, with an average value of 80 per pixel. For more in- 
formation about the stellar t argets, such as photom etry 
and radial velocity, please see lJacobson et aLl()2011l ) and 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

8 See http://www.astro.yale.edu/dokkum/lacosmic/ 
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Table 2 

Stellar Atmospheric Parameters & Fe abundances 



Cluster 


Star 


Teff 


log g 


vt 


[Fe I/H] 




# lines 


[Fe II/H] 


cr 


# Hn 






(K) 


(dex) 


(km s- 1 ) 


dex 


dex 




dex 


(dex) 




Be 18 


1163 


4650 


2.3 


1.15 


-0.30 


0.08 


67 


-0.31 


0.12 


8 


Be 18 


1383 


4400 


1.9 


1.30 


-0.34 


0.09 


71 


-0.33 


0.18 


8 


Be 21 


50 


4650 


2.1 


1.20 


-0.14 


0.09 


67 


-0.22 


0.08 


9 


Be 21 


51 


4520 


1.9 


1.40 


-0.28 


0.13 


70 


-0.26 


0.11 


9 


Be 22 


414 


4400 


2.1 


1.25 


-0.24 


0.12 


68 


-0.24 


0.12 


9 


Be 22 


6 13 


4000 


1.1 


1.50 


-0.30 


0.33 


61 


-0.26 


0.26 


6 


Be 17 


265 


4500 


2.1 


1.40 


-0.13 


0.14 


72 


-0.10 


0.16 


10 


Be 17 


569 


4300 


1.5 


1.50 


-0.11 


0.16 


73 


-0.11 


0.16 


11 


Be 17 


1035 


4300 


1.7 


1.50 


-0.13 


0.15 


73 


-0.19 


0.15 


11 


Be 32 


2 


4100 


1.0 


1.50 


-0.29 


0.14 


72 


-0.35 


0.27 


11 


Be 32 


4 


4100 


1.0 


1.50 


-0.33 


0.14 


72 


-0.41 


0.26 


11 


Be 32 


16 


4900 


2.7 


1.00 


-0.26 


0.10 


63 


-0.21 


0.15 


9 


Be 32 


18 


5000 


2.7 


1.30 


-0.19 


0.13 


66 


-0.21 


0.12 


9 


Be 39 


3 


4200 


1.6 


1.50 


-0.15 


0.12 


64 


—0.14 


0.14 


11 


Be 39 


5 


4450 


1.9 


1.40 


-0.11 


0.14 


69 


-0.16 


0.18 


11 


Be 39 


12 


4750 


2.3 


1.40 


-0.12 


0.17 


68 


-0.17 


0.07 


11 


Be 39 


1 1 


4750 


2.3 


1.40 


-0.15 


0.15 


70 


-0.13 


0.21 


11 


M 67 


105 


4450 


2.1 


1.50 


+0.03 


0.12 


64 


-0.01 


0.21 


11 


M 67 


141 


4700 


2.4 


1.50 


+0.09 


0.12 


65 


+0.07 


0.22 


11 


M 67 


170 


4300 


1.7 


1.50 


+0.02 


0.14 


60 


+0.00 


0.28 


11 


N 1193 


282 


4700 


2.2 


1.50 


-0.17 


0.14 


65 


-0.24 


0.16 


10 


N 1245 


10 


5000 


2.2 


1.20 


-0.01 


0.16 


69 


-0.02 


0.21 


11 


N 1245 


125 


4800 


2.5 


1.50 


+0.06 


0.18 


66 


+0.04 


0.12 


11 


N 1245 


160 


4800 


2.5 


1.30 


+0.03 


0.17 


63 


+0.03 


0.21 


11 


N 1245 


382 


4800 


2.3 


1.50 


+0.00 


0.13 


71 


-0.02 


0.15 


11 


N 1817 


73 


4850 


2.5 


1.40 


-0.04 


0.12 


73 


-0.07 


0.13 


11 


N 1817 


79 


5100 


2.6 


1.50 


-0.06 


0.14 


73 


-0.13 


0.11 


11 


N 188 


532 


4800 


2.9 


1.50 


+0.06 


0.19 


70 


-0.01 


0.11 


11 


N 188 


747 


4600 


2.6 


1.50 


+0.13 


0.19 


70 


+0.13 


0.21 


11 


N 188 


919 


4400 


2.2 


1.50 


+0.14 


0.19 


61 


+0.15 


0.16 


10 


N 188 


1224 


4700 


2.9 


1.50 


+0.13 


0.18 


71 


+0.12 


0.20 


11 


N 1883 


8 


4500 


1.5 


1.60 


-0.03 


0.13 


59 


-0.04 


0.20 


11 


N 1883 


9 


4600 


1.7 


1.60 


-0.04 


0.13 


61 


-0.06 


0.19 


11 


N 2141 


1007 


4100 


1.2 


1.60 


-0.09 


0.16 


59 


-0.07 


0.16 


11 


N 2141 


1348 


4100 


1.2 


1.50 


-0.09 


0.18 


47 


-0.08 


0.18 


8 


N 2158 


4230 


4400 


1.5 


1.50 


-0.05 


0.15 


58 


-0.11 


0.18 


11 


N 2194 


55 


5100 


2.5 


1.40 


-0.06 


0.14 


70 


-0.12 


0.17 


11 


N 2194 


57 


5100 


2.2 


1.40 


-0.06 


0.16 


70 


-0.11 


0.18 


11 


N 2355 


144 


5300 


3.3 


1.30 


-0.14 


0.21 


69 


-0.18 


0.22 


11 


N 2355 


398 


5200 


2.1 


1.90 


+0.05 


0.13 


71 


+0.01 


0.20 


11 


N 2355 


668 


5100 


2.6 


1.30 


-0.04 


0.13 


71 


-0.05 


0.15 


11 


N 6939 


121 


4200 


1.5 


1.60 


+0.00 


0.22 


61 


+0.02 


0.18 


11 


N 6939 


190 


5000 


2.6 


1.40 


+0.15 


0.17 


71 


+0.11 


0.21 


11 


N 6939 


212 


4300 


1.7 


1.50 


+0.03 


0.15 


65 


+0.07 


0.15 


11 


N 6939 


31 


4000 


0.9 


1.50 


-0.02 


0.18 


61 


+0.04 


0.26 


11 


N 7142 


196 


4500 


2.2 


1.40 


+0.08 


0.15 


55 


+0.10 


0.14 


11 


N 7142 


229 


4300 


1.7 


1.50 


+0.05 


0.15 


50 


+0.01 


0.12 


11 


N 7142 


377 


4450 


2.0 


1.50 


+0.08 


0.13 


56 


+0.03 


0.13 


11 


N 7142 


421 


4800 


2.5 


1.50 


+0.10 


0.14 


60 


+0.08 


0.14 


11 


PWM 4 


RGB 1 


4000 


1.2 


1.50 


-0.18 


0.16 


61 


-0.12 


0.30 


8 



IJacobson et all ([2007n (for NGC 6939). 

We carried out a detailed chemical abundance analy- 
sis of these stars using the same methods and line list 
as described here and in our previous work (see FJP10; 
updated log gf's for use with MOOG2010, as described 
below) for elements Fe, Na, Mg, Al, Si, Ca, Ti, Cr, Co, 
Ni, Zr, Ba, La and Eu. Here, we present the results for 
Fe and the neutron capture elements, and stellar atmo- 
spheric parameters; the other element abundances will be 
presented in a future paper (Jacobson et al., in prep.). 

All told, the total sample contains 50 stars in 19 clus- 
ters, with a maximum of four stars per cluster. Figure 
Q] shows portions of stellar spectra with absorption fea- 
tures of interest identified. Unfortunately, such small 
numbers of stars per cluster are insufficient to assess the 
true dispersions of element abundances in open clusters, 



and any conclusions we draw from their analysis are vul- 
nerable to the presence of "pathological" stars that are 
either non cluster members or just outliers. We have 
WIYN-Hydra single-order multi-object spectroscopy of 
stars in nine clusters in this study, as well as five others 
(| Jacobson et all 120071 1201 ll Friel et al. in prep.), that 
include La and Zr absorption features for as many as 30 
confirmed members per cluster. We will explore internal 
cluster abundance dispersions for La and Zr in a future 
paper. 



3. ANALYSIS 
3.1. Atmospheric parameters and [Fe/H] 
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Figure 2. Example syntheses of the La II 6262 feature in star PWM 4 RGB 1 (top panel) and the Eu II 6645 feature in star M 67 170 
(bottom panel). 



In thi s work , we used the 2010 version of MOOG 
(jSnedenl I1973D and spherical MARC^I model atmo- 
spheres (fGustafsso n et al.l [200 8V This is a change from 
our previous work where we exclusively used the 2002 
version of MOOG and an older grid of plane parallel 
MARCS models (|Bell et al.l [l976h . We note here that 
the change in the version of MOOG used in our analy- 
sis necessitated an update of the log gf's in our general 
line list (e.g., that used in FJP10), because they were de- 
termined differentially relative to Arcturus using MOOG 
and were therefore version-dependent. We have modified 
the log gf's of all the lines in our general list to be cor- 
rect for use with the 2010 version of MOOG and redeter- 
mined atmospheric parameters and element abundances 
for all the stars in our sample using the standard classi- 
cal abundance techniques (see, e.g., FJP10). Changes to 
log gf's were small with a mean difference of 0.01±0.03 
(1 cr) dex and a maximum of 0.11 dex for some Ti I 
lines. Relevant to this current work are the stellar atmo- 
spheric parameters, Fe and Zr abundances (see below). 
Updated detailed abundances of other elements for these 
stars are available by request to the authors. Differences 
in atmospheric parameters and element abundances were 
generally small - typically within 50 K, 0.1 dex, and 0.1 
km s _1 in T c ff, log g and v t , respectively (150 K, 0.3 



9 See htt p:/ /marcs. astro. uu.se/] the models were obtained from 
the website 2011 July. 



dex and 0.2 km s" 1 for the most extreme cases). These 
differences are well within the estimated uncertainties 
of the original analysis (see also Section 3.3). Tabled] 
presents the updated atmospheric parameters and iron 
abundances for the stars in this work. Note that [Fe/H] 
values are relative to a solar log e(Fe) = 7 .52, the default 
used in MOOG fsee ISneden et all 119911 ). and that "cr" 
denotes the (1 a) standard deviations of the Fe I and Fe 
II abundances. 

3.2. The Neutron capture elements 

We have employed both spectrum synthesis and equiv- 
alent width (EW) measurement of absorption features for 
the neutron-capture elements in this study. For Zr and 
Ba, we used EW measurements exclusively, while for La 
and Eu we used both EW measures and spectrum syn- 
thesis. Details of the analysis for each species are given 
below; here, however, we briefly summarize our methods 
of EW and spectrum synthesis analysis. 

Measurement of line EWs is straightforward and was 
carried out in the same way as done for Fe and the other 
elements in our previous work. That is, we fit Gaus- 
sians to the line profiles interactively in IRAF (though we 
also measured Ba features with Voigt profiles; see below). 
MOOG (version 2010) was used in our spectrum synthe- 
sis analysis, in a method identical to that employed to 
determine oxygen abundances in our previous work. We 
carefully matched the resolution of the synthetic spectra 
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Table 3 

Line List 



Element 


Wavelength 


E.P. 


log ftf 


log <H X ) 


i„„ ^ / "v" \ b 

log <H X ) 


i „ , /"v"^ c 
log e(X) 


[X/lcJ 




(A) 






Sun 


Sun 


Arcturus 


Arcturus 


Zr I 


6127.48 


0.150 


-1.180 


2.60 


2.94 


2.10 


+0.00 


Zr I 


6134.55 


0.000 


-1.426 


2.60 


3.02 


2.10 


+0.00 


Zr I 


6143.20 


0.070 


-1.252 


2.60 


2.90 


2.10 


+0.00 


Ba II 


5853.70 


0.604 


-1.010 s 


2.13 


2.20 


1.63 


+0.00 


Ba II 


6141.70 


0.704 


-0.070 e 


2.13 


2.59 


2.30 


+0.67 


Ba II 


6496.90 


0.604 


-0.377° 


2.13 


2.41 


2.06 


+0.43 


La II 


6262.29 


0.403 


-1.220 6 


1.13 


1.22 


0.45 


-0.18 


La II 


6390.48 


0.320 


-1.410 c 


1.13 


1.20 


0.49 


-0.14 


Eu II 


6437.64 


1.319 


-0.320 6 


0.52 


0.53 


0.25 f 


+0.23 


Eu II 


6645.10 


1.379 


+0.120 c 


0.52 


0.49 


0.25 


+0.23 



a Literature values as described in text. 

b Values found via EW measures (Zr, Ba) or spectrum synthesis (La, Eu) in this work, as 
described in text. 
c This study. 

d Arcturus [X/Fe] ratios calculated relative to literature solar values, adopting [Fe/H] = —0.50. 
e The total log gf for the feature ; log gf's for the individu al components in La and Eu are 
presented in ILaw lcr et al. (2001a|) and ILaw ler ct al. (2001b), respe ctivel y. T otal log gf's and 
those of individual components for Ba are adopted from McWilliam ( 1998) and Prochaska ct al. 
lt2000l V 

f This abundance is forced to match that for the Eu II A6645 feature; see details in the text. 



to the observed spectra by synthesizing absorption fea- 
tures of other elements near the feature of interest whose 
abundances we determined via EW measurement (e.g., 
Fe, Ca, or Ni lines). Next, we generated sets of three syn- 
thetic spectra at a time with the abundance of the ele- 
ment of interest varied in a stepsize 0.30 dex and the best 
fit abundance was determined by eye. The uncertainty 
in the best fit abundance was then found by decreasing 
the step size between the three synthetic spectra, and 
then varying the level of the continuum by its approxi- 
mate uncertainty, until no single synthetic spectra could 
be identified as the best. 

3.2.1. Zirconium 

The zirconium abundances are based on equivalent 
width (EW) measurements of three Zr I lines at A6127, 
6134, and 6143 A. The log gf's of these li nes were de- 
termin ed by measuring their EWs in the iHinkle et al.l 
( 2000) spectrum of Arcturus a nd forcing the a bundances 
to match the Zr abundance of IPeterson et al.l (1993; log 
e(Zr) = 2.10, normalized to the solar abundance value of 
lAnders &: Grevesselll989L 2.60). The atomic information 
for these lines are given in Table 0B 

Our log gf values for these Zr lines are generally smaller 
than those found in t he literature by 0.05 to 0.15 dex 
(e.g. JYong et al.l 120051: Dohnson et al.ll2012r i. As a check 
on what solar abundances we would derive using these 
transition probabilities, we measured these lines i n the 
solar flux atlas that accompanies the IHinkle et al.l Arc- 
turus atlas. The solar abundances we derived are shown 
in Table [3l and are ^0.30-0.40 dex larger than litera- 
ture values (e.g..lAnders k, GrevesseTl 989; As plund et al.l 
[20091: ILodders et al.ll2009ir This indicates that our choice 
of log gf is possibly responsible for the systematically 

10 Technically, five isotopes of Zr contribute to its total 
abundance in solar syste m: 90 Zr, 91 Zr, 92 Zr, 94 Zr, and 96 Zr 
IjAnders fe Grevess c 1989j). Only odd isotopes are subject to hy- 
pcrfinc splitting, and as 91 Zr comprises 11.22% of the solar system 
abundance (and therefore only a fraction of the Zr line profiles), 
hyperfine and isotopic splitting are not considered here. 



higher Zr abundances compared to literature values for 
many stars in our sample, as described below, though we 
note that our log gf's for Zr 6134 and 6 143 agree within 
0.1 dex of those of Joh nson et al.l (|2012l ). who also found 
[Zr/Fe] = 0.00 for Arcturus, a fter a dopting log e(Zr) = 
2.60 for the Sun. lYong et all (|2005f) . whose log gf's for 
these Zr lines are 0.15 dex larger than ours, found log 
N(Zr) = 2.80 for the Sun and [Zr/Fe] = -0.27 for Arc- 
turus. iSmith et al.l (12000T). wh ose log gf's were identical 
to those of lYong et all (j2005h . found log N(Zr) = 2.88 
for the Sun, and [Zr/Fe] = -0.24 for Arcturus. Clearly, 
differences in log gf's do not fully account for differences 
in abundances here, which, in general, should be kept in 
mind when comparing results of different studies in the 
literature. 

Note that for this present work, the EWs of A6134 
and 6143A Zr I lines have been remeasured in all stars 
from our previous papers and therefore may be differ- 
ent from values published there. The abundances from 
the three Zr I lines are in general good agreement with 
abundance differences smaller than 0.1 dex for all but the 
lower S/N stars (e.g., NGC 2355 stars). For most of the 
stars from our previous work, remeasurement of A6134 
and 6143 features, along with the addition of A6127 line 
decreased the line-by-line dispersion in abundances, and 
there are no systematic differences between the individ- 
ual line abundances. However, for a handful of stars, 
the abundance from one line differed greatly from the 
other two; in these cases, visual inspection of the spec- 
tra and remeasurement of the EWs did not resolve the 
discrepency or identify its source. In such cases, we have 
included lines with discrepant abundances, even though 
the resulting dispersion in the mean abundance is larger. 
Individual star [Zr/Fe] abundances, calculated relative to 
the [Fe I/H] ratios given in Table [2] are shown in Table 
[4j Mean cluster [Zr/Fe] abundances, with (1 a) standard 
deviations, are shown in Table [H 

3.2.2. Barium 
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Table 4 

Individual Star [Zr/Fe] Ratios 



Table 5 

Individual Star [Ba/Fe] Ratios 







Zr A6127 


Zr A6134 


Zr A6143 


[Zr/Fe] 


[Zr/Fe] 








Ba A5853 


Ba A6141 


Ba A6496 


[Ba/Fe] 


[Ba/Fe] 


Cluster 


Star 


EW 


EW 


EW 


mean 




<T 


Cluster 


Star 


EW 


EW 


EW 


mean 


<7 


Be 17 


265 


-0.09 


-0.18 


-0.10 


— 


.12 


0. 


.05 


Be 


17 


265 


+0.04 


-0.24 


+0.05 


+0.04 


0.01 


Be 17 


569 


-0.15 


-0.42 


-0.23 


— 0. 


.27 


0. 


.14 


Be 


17 


569 


+0.09 


-0.23 


+0.09 


+0.09 


0.00 


Be 17 


1035 


+0.20 


-0.14 


-0.02 


+0. 


.01 





.17 


Be 


17 


1035 


+0.01 


-0.16 


-0.03 


-0.01 


0.03 


Be 18 


1163 


+0.20 


+0.04 


+0.05 


+0. 


.10 





.09 


Be 


18 


1163 


+0.42 


+0.11 


+0.36 


+0.39 


0.04 


Be 18 


1383 


+0.06 


+0.02 


+0.03 


+0 


.04 


0. 


.02 


Be 


18 


1383 


+0.47 


+0.16 


+0.37 


+0.42 


0.07 


Be 21 


50 


+0.24 


+0.17 


+0.17 


+0 


.19 





.0 1 


Be 


21 


50 


+0.64 


+0.20 


+0.49 


+0.56 


0.11 


Be 21 


51 


+0.10 


+0.01 


+0.06 


+0 


.06 





.05 


Be 


21 


51 


+0.45 


+0.26 


+0.40 


+0.43 


0.04 


Be 22 


414 


+0.27 


+0.06 


+0.10 


+0. 


.14 


0. 


11 


Be 


22 


414 


+0.51 


+0.08 


+0.38 


+0.45 


0.09 


Be 22 


643 


+0.80 


+0.71 


+0.94 


+0 


.82 





.12 


Be 


22 


643 


+1.01 


+0.41 


+0.51 


+0.76 


0.35 


Be 32 


2 


-0.08 


+0.05 


-0.01 


-0. 


.01 


0. 


,07 


Be 


32 


2 


+0.30 


-0.11 


+0.08 


+0.19 


0.16 


Be 32 


4 


-0.14 


-0.13 


-0.09 


-0. 


.12 





.03 


Be 


32 


4 


+0.28 


-0.14 


+0.09 


+0.18 


0.13 


Be 32 


16 


+0.07 


+0.21 


+0.18 


+0 


.15 





.07 


Be 


32 


16 


+0.37 


-0.17 


+0.16 


+0.27 


0.15 


Be 32 


18 


+0.21 


+0.23 


+0.28 


+0 


.24 





.04 


Be 


32 


18 


+0.32 


+0.01 


+0.12 


+0.22 


0.14 


Be 39 


3 


-0.05 


-0.15 


+0.08 


-0 


.04 


0. 


.12 


Be 


39 


3 


+0.24 


-0.19 


+0.08 


+0.16 


0.11 


Be 39 


5 


+0.07 


-0.18 


-0.03 


-0 


.05 





.13 


Be 


39 


5 


+0.03 


-0.55 


+0.04 


+0.04 


0.01 


Be 39 


12 


+0.05 


+0.09 


-0.05 


+0. 


.03 





.07 


Be 


39 


12 


+0.05 


-0.17 


+0.09 


+0.07 


0.03 


Be 39 


11 


+0.16 


+0.37 


+0.05 


+0 


.19 


0. 


.16 


Be 


39 


14 


-0.07 


-0.16 


+0.12 


+0.02 


0.13 


M 67 


105 


+0.10 


-0.00 


+0.10 


+0 


.07 





.06 


Al 


67 


105 


+0.07 


-0.31 


+0.00 


+0.04 


0.05 


M 67 


141 


-0.03 


-0.06 


-0.02 


-0 


.04 


0. 


.02 


Al 


67 


111 


+0.05 


-0.30 


+0.14 


+0.10 


0.06 


M 67 


170 


-0.06 


-0.11 


+0.01 


— 0. 


.05 


0. 


.06 


Al 


67 


170 


+0.19 


-0.17 


+0.14 


+0.17 


0.04 


N 1193 


282 


+0.48 


+0.06 


+0.45 


+0 


.33 





.23 


N 1193 


282 


+0.15 


-0.05 


+0.13 


+0.14 


0.01 


N 1245 


10 


+0.20 


+0.11 


+0.27 


+0. 


.19 





.08 


N 1245 


10 


+0.47 


-0.08 


+0.23 


+0.35 


0.17 


N 1245 


125 


-0.02 


+0.02 


+0.22 


+0 


.07 


0. 


.13 


N 1245 


125 


+0.70 


+0.15 


+0.31 


+0.51 


0.28 


N 1245 


160 


+0.54 


+0.55 


+0.52 


+0 


.54 





.02 


N 1245 


160 


+0.43 


+0.20 


+0.43 


+0.43 


0.00 


N 1245 


382 


+0.26 


+0.25 


+0.29 


+0. 


.27 





.02 


N 1245 


382 


+0.44 


+0.05 


+0.34 


+0.39 


0.07 


N 1817 


73 


+0.34 


+0.17 


+0.17 


+0 


.23 





.10 


N 1817 


73 


+0.56 


+0.15 


+0.15 


+0.36 


0.29 


N 1817 


79 


+0.27 


-0.02 


+0.11 


+0. 


.12 


0. 


.15 


N 1817 


79 


+0.32 


+0.04 


+0.04 


+0.18 


0.20 


N 188 


532 


+0.28 


+0.06 


+0.45 


+0 


.26 





.20 


N 188 


532 


-0.17 


-0.39 


-0.10 


-0.14 


0.05 


N 188 


717 


+0.03 


-0.07 


+0.04 


-0 


.00 





.06 


N 188 


747 


-0.12 


-0.34 


+0.14 


+0.01 


0.18 


N 188 


919 


+0.07 


-0.27 


-0.10 


-0 


.10 


0. 


.17 


N 188 


919 


-0.14 


-0.25 


+0.18 


+0.02 


0.23 


N 188 


1224 


+0.13 


+0.13 


+0.19 


+0 


.15 





.03 


N 188 


1224 


+0.33 


-0.07 


+0.22 


+0.28 


0.08 


N 1883 


8 


+0.05 


-0.28 


-0.15 


-0 


.13 





.17 


N 1883 


8 


+0.56 


+0.13 


+0.34 


+0.45 


0.16 


N 1883 


9 


-0.01 


-0.11 


-0.13 


-0. 


.08 


0. 


.06 


N 1883 


9 


+0.58 


+0.20 


+0.35 


+0.46 


0.16 


N 2141 


1007 


-0.08 


-0.04 


+0.02 


-0 


.03 





.05 


N 2141 


1007 


+0.37 


-0.01 


+0.17 


+0.27 


0.14 


N 2141 


1348 


+0.89 


+0.45 


+0.85 


+0 


.73 





.21 


N 2141 


1348 


+0.57 


+0.33 


+0.51 


+0.54 


0.04 


N 2158 


4230 


-0.03 


-0.10 


-0.02 


— 0. 


.05 


0. 


.01 


N 2158 


4230 


+0.36 


+0.04 


+0.18 


+0.27 


0.13 


N 2194 


55 


+0.15 


+0.15 


+0.29 


+0 


.20 





.08 


N 2194 


55 


+0.68 


+0.13 


+0.44 


+0.56 


0.17 


N 2194 


57 


+0.17 


-0.16 


+0.36 


+0. 


.12 


0. 


.26 


N 2194 


57 


+0.63 


+0.20 


+0.37 


+0.50 


0.18 


N 2355 


144 


+0.92 




+0.59 


+0 


.76 


0. 


.23 


N 2355 


144 


+0.47 


+0.17 


+0.46 


+0.46 


0.01 


N 2355 


398 


+0.19 


+0.36 


+0.60 


+0 


.38 





21 


N 2355 


398 


+0.90 


+0.44 


+0.75 


+0.82 


0.11 


N 2355 


668 


+0.23 


+0.47 


+0.28 


+0 


.33 


0. 


.13 


N 2355 


668 


+0.52 


+0.17 


+0.40 


+0.46 


0.08 


N 6939 


121 


-0.70 


-0.66 


-0.62 


-0 


.66 


0. 


.01 


N 6939 


121 


-0.16 


-0.30 


-0.34 


-0.25 


0.13 


N 6939 


190 


+0.55 


+0.48 


+0.82 


+0 


.62 





.18 


N 6939 


190 


+0.98 


+0.46 


+0.60 


+0.79 


0.27 


N 6939 


212 


+0.08 


-0.17 


-0.13 


-0 


.07 


0. 


.13 


N 6939 


212 


+0.38 


+0.01 


-0.05 


+0.17 


0.30 


N 6939 


31 


+0.28 


+0.11 


+0.21 


+0 


.20 


0. 


.09 


N 6939 


31 


+0.53 


+0.07 


+0.04 


+0.29 


0.35 


N 7142 


196 


+0.17 


-0.03 


+0.10 


+0 


.08 


0. 


.10 


N 7142 


196 


+0.19 


-0.23 


+0.18 


+0.18 


0.01 


N 7142 


229 


+0.29 


+0.04 


+0.20 


+0 


.18 


0. 


.13 


N 7142 


229 


+0.24 


-0.09 


+0.21 


+0.22 


0.02 


N 7142 


377 


+0.17 


+0.05 


+0.21 


+0 


.14 





.08 


N 7142 


377 


+0.21 


-0.18 


+0.19 


+0.20 


0.01 


N 7142 


421 


+0.14 


-0.10 


+0.24 


+0 


.09 





.17 


N 7142 


421 


+0.20 


-0.22 


+0.16 


+0.18 


0.03 


PWM 4 


RGB 1 


+0.23 


+0.11 


+0.22 


+0 


.19 





.07 


PWM 4 


RGB 1 


+0.44 


+0.08 


+0.24 


+0.34 


0.14 



Three barium lines were measured for determination 
of abundances: A5853, 6141, and 6496 A. All three lines 
are quite strong and broad, and the wings of some are 
blended with other spectral features (namely A6496). 
Barium features are subject to both hyperfine and iso- 
topic broadening from seven different isotopes: 130 Ba, 
132 Ba, 134 Ba, 135 Ba, 13 6 Ba, 1 37 Ba, and 138 Ba. According 
to lAnders k. Grevessd (|1989D . these isotopes contribute 
0.106%, 0.101%, 2.417%, 6.592%, 7.854%, 11.23%, and 
71.70% of the total solar system Ba abunda nce, respec- 
tively. We hav e employed the l i ne lis ts of iMcWilliaml 
(1998) (see also iProchaska et al.l 120001 ) . which incorpo- 
rate these isotopic fractions within the log gf's, as well 
as the s-process and r-process contributions to the solar 
system barium abundance (85% and 15%, respectively). 

Abundances from these lines were calculated using 



Note. — Individual line [Ba/Fc] ratios are calculated differentially relative to 
the solar Ba abundances listed in Tablc[3] The mean [Ba/Fc] values are calculated 
using the 5853 and 6496 features only. See the text for more information. 

measured EWs. To determine the best me thod to mea- 
sure E Ws, we measured these lines in the IHinkle et al.l 
(2000) Arcturus and solar flux atlases by fitting both 
Gaussian and Voigt profiles to them. The Ba abun- 
dances corresponding to the Gaussian line profiles are 
indicated in Table [3] Voigt profiles resulted in abun- 
dances 0.2-0.6 dex larger than that those for the Gaus- 
sian fits, and much more deviant f rom the solar value 
found in the literature (e.g. , 2.13 - lAnders fc Grevessel 
119891 2.24 - iGrevesse et al.l I1996T ). We also measured 
Ba abundances in the Sun from these lines via spectrum 
synthesis and find values more consistent with the solar 
value: log e(Ba) = 2.13, 2.23, and 2.26 for A5853, 6141, 
and 6496, respectively. Using the same three Ba lines 
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Table 6 

Individual Star [La/Fe] Ratios 



Cluster 


otar 


[±ja/ _rC| AOZOZ 


unc. AOZOZ 


1 Liu./ r cj AOoyu 


unc. ADoyu 


mean [La/ Pe] 


_[T „ IT?. 


Be 17 


265 


-0.33 


0.30 


-0.03 


0.10 


-0.18 


0.21 


Be 17 


569 


-0.30 


0.10 


-0.30 


0.10 


-0.30 


0.00 


Be 17 


1035 


-0.23 


0.07 


-0.08 


0.10 


-0.15 


0.11 


Be 18 


1163 


+0.16 


0.05 


+0.17 


0.05 


+0.17 


0.01 


Be 18 


1383 


+0.03 


0.10 


+0.18 


0.06 


+0.10 


0.11 


Be 21 


50 


+0.10 


0.05 


+0.13 


0.06 


+0.12 


0.02 


Be 21 


51 


+0.14 


0.05 


+0.17 


0.05 


+0.16 


0.02 


Be 22 


414 


+0.16 


0.05 


+0.20 


0.05 


+0.18 


0.03 


Be 22 


643 


+0.07 


0.05 


+0.36 


0.07 


+0.21 


0.21 


Be 32 


2 


-0.24 


0.05 


-0.12 


0.06 


-0.18 


0.08 


Be 32 


4 


-0.33 


0.10 


-0.18 


0.07 


-0.26 


0.11 


Be 32 


16 


+0.10 


0.06 


+0.12 


0.07 


+0.11 


0.01 


Be 32 


18 


+0.01 


0.06 


+0.03 


0.06 


+0.02 


0.01 


Be 39 


3 


-0.23 


0.06 


-0.08 


0.06 


-0.15 


0.11 


Be 39 


5 


-0.20 


0.10 


+0.02 


0.06 


-0.09 


0.16 


Be 39 


12 














Be 39 


14 


-0.06 


0.10 


+0.04 


0.08 


-0.01 


0.07 


M 67 


105 


-0.14 


0.06 


-0.09 


0.05 


-0.12 


0.04 


M 67 


111 


-0.15 


0.10 


-0.13 


0.07 


-0.14 


0.01 


M 67 


170 


-0.25 


0.08 


-0.15 


0.09 


-0.20 


0.07 


N 1193 


282 


+0.01 


0.08 


+0.11 


0.07 


+0.06 


0.07 


N 1245 


10 


-0.35 


0.15 


-0.30 


0.10 


-0.32 


0.04 


N 1245 


125 


+0.03 


0.10 


-0.06 


0.10 


-0.01 


0.06 


N 1245 


160 


+0.21 


0.15 




0.00 


+0.21 


0.15 


N 1245 


382 


+0.09 


0.10 


+0.09 


0.12 


+0.09 


0.00 


N 1817 


73 


+0.10 


0.10 


+0.13 


0.08 


+0.12 


0.02 


N 1817 


79 


+0.00 


0.10 


+0.15 


0.10 


+0.07 


0.11 


N 188 


532 


+0.03 


0.10 


-0.10 


0.10 


-0.04 


0.09 


N 188 


747 


-0.19 


0.10 


-0.12 


0.10 


-0.15 


0.05 


N 188 


919 


-0.33 


0.10 


-0.23 


0.10 


-0.28 


0.07 


N 188 


1224 


+0.18 


0.08 


+0.30 


0.12 


+0.24 


0.08 


N 1883 


8 


-0.15 


0.08 


-0.10 


0.16 


-0.12 


0.04 


N 1883 


9 


+0.01 


0.07 


-0.01 


0.10 


+0.00 


0.01 


N 2141 


1007 


-0.17 


0.05 


-0.06 


0.10 


-0.12 


0.08 


N 2141 


1348 


+0.08 


0.04 


+0.18 


0.06 


+0.13 


0.07 


N 2158 


4230 


-0.16 


0.07 


-0.15 


0.08 


-0.15 


0.01 


N 2194 


55 


+0.13 


0.10 


+0.06 


0.10 


+0.10 


0.05 


N 2194 


57 


+0.00 


0.08 


-0.05 


0.10 


-0.03 


0.04 


N 2355 


144 


+0.28 


0.10 


+0.38 


0.10 


+0.33 


0.07 


N 2355 


398 


+0.07 


0.08 


+0.10 


0.08 


+0.09 


0.02 


N 2355 


668 


+0.10 


0.10 


+0.13 


0.10 


+0.12 


0.02 


N 6939 


121 


-0.41 


0.10 


-0.39 


0.10 


-0.40 


0.01 


N 6939 


190 


+0.34 


0.10 


+0.49 


0.10 


+0.42 


0.11 


N 6939 


212 


-0.19 


0.08 


-0.12 


0.09 


-0.15 


0.05 


N 6939 


31 


-0.39 


0.05 


-0.29 


0.12 


-0.34 


0.07 


N 7142 


196 


-0.24 


0.10 


-0.19 


0.10 


-0.21 


0.04 


N 7142 


229 


-0.19 


0.07 


-0.06 


0.07 


-0.12 


0.09 


N 7142 


377 


-0.21 


0.07 


-0.09 


0.06 


-0.15 


0.08 


N 7142 


421 


-0.16 


0.06 


-0.16 


0.08 


-0.16 


0.00 


PWM 4 


RGB 1 


-0.01 


0.05 


+0.10 


0.06 


+0.05 


0.08 



as this study, [Smith et al l (12000) found [Ba/Fe] = -0.05 
for Arcturus; lYong et al.l (|2005D found [Ba/Fe] = +0.09 
based on the A5853 line only. Our value for A5853 is 
consistent with these values. 

As a result of this analysis of the Sun and Arcturus, 
we measured all three Ba features in our sample stars 
by fitting Gaussian profiles, as these values are more 
consistent with those both from the spectrum synthe- 
sis analysis and from the literature. Given the especially 
large discrepancy between the spectrum synthesis and 
EW abundance for the A6141 line in the Sun and Arc- 
turus, and that its abundance is also often discrepant 
from those of the other Ba lines in our program stars, we 
have opted not to include this line in the calculation of 
mean Ba abundances for our stellar sample, though we 
include it in our results for completeness. Table [5] shows 
the [Ba/Fe] ratios for each star calculated from the in- 



dividual lines, with the mean and standard deviations 
calculated from the A5853 and 6496 lines only. Note that 
these [Ba/Fe] ratios are calculated relative to our deter- 
mined solar Ba abundances in column 6 of Table Q to 
minimize systematics due to the zero-point of our abun- 
dance scale. Mean cluster [Ba/Fe] abundances, with (1 
a) standard deviations, are shown in Table [8] 

3.2.3. Lanthanum 

139 L a comprises 99.911% of th e solar system La abun- 
dance (|Anders fc G rcvessc 1989), and so La absorption 
features are subject to hyperfine broadening . We have 
adopted the line list of lLawler et al.l (|2001afl , which in- 
cludes data for individual hyperfine lines, for the analysis 
of the A6262 and A6390 A features. We have determined 
La abundances for all stars in this study via both spec- 
trum synthesis and EW measurement of the lines. We 
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Table 7 

Individual Star [Eu/Fe] Ratios 



Cluster 


Diar 


l"FTn /TTnl \&A17 
[H/U j r CJ AD4o I 


uric. AD^to ( 


[HjU/ rej ado^q 


unc. A6645 


nica,n [Eu/Fc] 




Be 17 


265 


+0.02 


0.12 


+0.02 


0.08 


+0.02 


0.00 


Be 17 


569 


-0.15 


0.15 


-0.15 


0.08 


-0.15 


0.00 


Be 17 


1035 


+0.02 


0.13 


-0.03 


0.07 


-0.00 


0.04 


Be 18 


1163 


+0.11 


0.15 


+0.17 


0.07 


+0.14 


0.04 


Be 18 


1383 




0.00 


+0.16 


0.07 


+0.16 


0.07 


Be 21 


50 


-0.07 


0.15 


-0.12 


0.10 


-0.10 


0.04 


Be 21 


51 


+0.14 


0.08 


+0.05 


0.08 


+0.10 


0.06 


Be 22 


414 


+0.23 


0.13 


+0.18 


0.08 


+0.21 


0.04 


Be 22 


643 


+0.41 


0.08 


+0.24 


0.08 


+0.32 


0.12 


Be 32 


2 


-0.14 


0.10 


-0.02 


0.06 


-0.08 


0.08 


Be 32 


4 


+0.00 


0.08 


+0.04 


0.06 


+0.02 


0.03 


Be 32 


16 


+0.05 


0.10 


+0.08 


0.08 


+0.07 


0.02 


Be 32 


IS 


+0.13 


0.07 


+0.08 


0.05 


+0.11 


0.04 


Be 39 


3 


-0.26 


0.12 


-0.09 


0.07 


-0.17 


0.12 


Be 39 


5 


+0.00 


0.13 


-0.13 


0.07 


-0.07 


0.09 


Be 39 


12 


+0.01 


0.10 


-0.09 


0.15 


-0.04 


0.07 


Be 39 


14 


-0.01 


0.13 


-0.01 


0.08 


-0.01 


0.00 


M 67 


105 


-0.14 


0.10 


-0.07 


0.07 


-0.11 


0.05 


M 67 


141 


-0.10 


0.10 


-0.05 


0.07 


-0.08 


0.04 


M 67 


170 


-0.28 


0.10 


-0.13 


0.07 


-0.21 


0.11 


N 1193 


282 


+0.06 


0.10 


-0.04 


0.10 


+0.01 


0.07 


N 1245 


10 














N 1245 


125 














N 1245 


160 














N 1245 


382 














N 1817 


73 


-0.12 


0.15 


-0.07 


0.10 


-0.10 


0.04 


N 1817 


79 


-0.20 


0.15 


-0.10 


0.10 


-0.15 


0.07 


N 188 


532 


+0.08 


0.15 


-0.22 


0.10 


-0.07 


0.21 


N 188 


747 


-0.44 


0.20 


-0.04 


0.12 


-0.24 


0.28 


N 188 


919 


-0.53 


0.20 


-0.13 


0.10 


-0.33 


0.28 


N 188 


1224 


-0.34 


0.20 


+0.16 


0.10 


-0.09 


0.35 


N 1883 


8 


-0.30 


0.13 


-0.30 


0.15 


-0.30 


0.00 


N 1883 


9 


-0.09 


0.10 


-0.20 


0.10 


-0.15 


0.08 


N 2141 


1007 


-0.22 


0.10 


-0.14 


0.08 


-0.18 


0.06 


N 2141 


1348 


-0.22 


0.10 


-0.12 


0.10 


-0.17 


0.07 


N 2158 


4230 


-0.16 


0.10 


-0.11 


0.10 


-0.14 


0.04 


N 2194 


55 


-0.30 


0.18 


-0.10 


0.10 


-0.20 


0.14 


N 2194 


57 


-0.30 


0.15 


-0.30 


0.15 


-0.30 


0.00 


N 2355 


144 


+0.28 


0.15 




0.00 


+0.28 


0.15 


N 2355 


398 


-0.01 


0.10 


-0.16 


0.10 


-0.09 


0.11 


N 2355 


668 


-0.35 


0.30 


-0.17 


0.10 


-0.26 


0.13 


N 6939 


121 


-0.31 


0.15 


-0.21 


0.10 


-0.26 


0.07 


N 6939 


190 






-0.13 


0.10 


-0.13 


0.10 


N 6939 


212 


-0.34 


0.17 


-0.11 


0.08 


-0.23 


0.16 


N 6939 


31 


-0.26 


0.15 


-0.34 


0.10 


-0.30 


0.06 


N 7142 


196 


-0.29 


0.10 


-0.19 


0.10 


-0.24 


0.07 


N 7142 


229 


-0.46 


0.15 


-0.21 


0.10 


-0.34 


0.18 


N 7142 


377 


-0.54 


0.30 


-0.19 


0.10 


-0.36 


0.25 


N 7142 


421 


-0.31 


0.15 


-0.21 


0.10 


-0.26 


0.07 


PWM 4 


RGB 1 


+0.07 


0.10 


+0.07 


0.08 


+0.07 


0.00 



take the spectrum synthesis abundance values to be more 
reliable, since fitting the line features with Gaussians to 
determine an equivalent width may not necessarily well 
fit the broad La features in some stars. Figure [2] shows 
an example spectrum synthesis of the La II 6262 feature 
in one of our program stars. 

First though, it is useful to confirm the zero-point 
of our La abundances. To do so, we have determined 
abundances of the Sun a nd Arcturus v i a spe ctrum syn- 
thesis, again using the iHinkle et al.l ([2000) Arcturus 
atla s and the a c compan ying solar flux atlas. Using 
the lLawler et al.l (|2001aD line list, we find log e(La) = 
1.22 and 1.20 for A6262 and A6390 A for the Sun, re- 
spectively (see Table EJ). These values, while slightly- 
larger t han the log e(La) = 1.13 found bv lLawler et al.l 
(l2001alK are in good agreem ent with solar values found by 
lAnders &; Grevessi (|1989f) and meteoritic values found 



bv lLodders et a l. ( 2009]E3- The individual line La abun- 
dances for Arcturus also agree very well (within 0.04 dex 
of each other); however, our mean [La/Fe] ratio for Arc- 
tur us, — 0.16±0.0 3, is .10 dex lower than that found 
by Uohnson et all (|2012t ) using the same line list and 
spectrum. Our v alue is even lower than that found by 
lYong et alJ (|2005D and lSmith et al.l (|2000H for Arcturus: 
[La/Fe] = +0.03±0.09 and [La/Fe] = +0.21, respectively 
(though we note that this difference in [La/Fe] values is 
largely driven by a ^0.2 dex difference in [Fe/H] found 
for Arcturus between lSmith et al.l (12000T) and other stud- 
ies cited here). We adopt the lLawler et al . (2001a) solar 
La abundance of 1.13 to calculate [La/Fe] ratios for this 

11 As pointed out by the referee, Lawlcr ct al. (2001a) used the 
empir ical model atmosphere of the Sun from Holwcg er He Muellerl 
H1974T ) in their analysis. This may at least partly explain the dif- 
ferences between our determined La abundances and their values. 
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Table 8 

Cluster mean abundances 



Cluster 


[Fe/H] 


(7 


[Zr/Fe] 


(7 


[Ba/Fe] 


(7 


[La/Fe] 


a 


[Eu/Fe] 


<7 


# St; 


Be 17 


-0.12 


0.01 


-0.13 


0.14 


+0.04 


0.05 


-0.21 


0.08 


-0.04 


0.09 


3 


Be 18 


-0.32 


0.03 


+0.07 


0.04 


+0.41 


0.02 


+0.14 


0.05 


+0.15 


0.01 


2 


Be 21 


-0.21 


0.10 


+0.13 


0.09 


+0.50 


0.09 


+0.14 


0.03 


+0.00 


0.14 


2 


Be 22 1 


-0.24 


0.12 


+0.14 


0.11 


+0.45 


0.09 


+0.18 


0.03 


+0.21 


0.04 


1 


Be 32 


-0.27 


0.06 


+0.07 


0.16 


+0.22 


0.04 


-0.08 


0.17 


+0.03 


0.08 


4 


Be 39 


-0.13 


0.02 


+0.03 


0.11 


+0.07 


0.06 


-0.08 


0.07 


-0.07 


0.07 


4 


M 67 


+0.05 


0.04 


-0.01 


0.07 


+0.10 


0.07 


-0.15 


0.04 


-0.13 


0.07 


3 


NGC 188 


+0.12 


0.04 


+0.08 


0.16 


+0.04 


0.17 


-0.06 


0.22 


-0.18 


0.12 


1 


NGC 1193 


-0.17 


0.13 


+0.33 


0.23 


+0.14 


0.01 


+0.06 


0.07 


+0.01 


0.07 


1 


NGC 1245 


+0.02 


0.03 


+0.27 


0.20 


+0.42 


0.07 


-0.01 


0.23 






4 


NGC 1817 


-0.05 


0.01 


+0.18 


0.08 


+0.27 


0.13 


+0.10 


0.04 


-0.13 


0.04 


2 


NGC 1883 


-0.04 


0.02 


-0.11 


0.04 


+0.46 


0.02 


-0.06 


0.08 


-0.23 


0.11 


2 


NGC 2141 


-0.09 


0.01 


+0.50 


0.71 


+0.41 


0.19 


+0.01 


0.18 


-0.18 


0.01 


2 


NGC 2158 


-0.05 


0.14 


-0.05 


0.04 


+0.27 


0.13 


-0.15 


0.01 


-0.14 


0.04 


1 


NGC 2194 


-0.06 


0.01 


+0.16 


0.06 


+0.41 


0.19 


+0.04 


0.09 


-0.25 


0.07 


2 


NGC 2355 


-0.04 


0.10 


+0.49 


0.24 


+0.58 


0.21 


+0.18 


0.13 


-0.02 


0.28 


3 


NGC 6939 2 


+0.00 


0.03 


-0.31 


0.31 


+0.07 


0.28 


-0.30 


0.13 


-0.26 


0.04 


3 


NGC 7142 


+0.08 


0.02 


+0.12 


0.05 


+0.20 


0.02 


-0.16 


0.04 


-0.30 


0.06 


1 


PWM 4 


-0.18 


0.16 


+0.19 


0.07 


+0.34 


0.14 


+0.05 


0.08 


+0.07 


0.01 


1 



1 [X/Fe] ratios calculated without star 643. 

2 [X/Fe] ratios calculated without star 190. 

study. Individual star [La/Fe] ratios as determined via 
spectrum synthesis are shown in Table[6] The abundance 
uncertainty for each line given in this Table represents 
the measurement uncertainty in the synthesis technique, 
as described in the introduction of this Section. Mean 
cluster [La/Fe] abundances, with (1 a) standard devi- 
ations, are shown in Table |S] (For cluster abundances 
based on the measures of a single star only, the 1 a 
value represents the uncertainty (standard deviation) in 
its abundance from Table 

Abundances from each La II feature were also deter- 
mined using measured EWs and the blends driver in 
MOOG, to correctly account for the hyperfine structure 
of the lines. Differences between synthesis and EW abun- 
dances for each line were plotted against [Fe/H], T e g, 
and log g to look for systematic trends. For the 6262 
feature, the EW abundances were ~0.07 smaller than 
the synthesis abundances, with a standard deviation of 
^0.15 dex. While no trend with [Fe/H] is present, the 
difference between EW and synthesis abundances grows 
larger (that is, EW abundances get increasingly smaller 
relative to the synthesis values) with increasing T e g and 
log g. For the 6390 feature, the difference between EW 
and synthesis abundances (in the sense EW — synthesis) 
is ~ +0.13, with a standard deviation of ^0.12 dex. No 
trend is present versus [Fe/H], T e ff, or log g. 

3.2.4. Europium 

Like lanthanum, europium is subject to strong hyper- 
fine splitting. However, while La features are predom- 
inantly composed of a single isotope, Eu features are 
also subject to splitting between two isotopes. We have 
adopted the l ine list for the Eu II A6437 and A6645 A 
features from lLawler et al.l (|2001bf ) , assuming the stan- 
dard s olar system isotopic mix of 47.8% 151 Eu and 52.2% 
153 Eu (jAnders fc GrevessdH989f ). 

It is important to state that, according to lLawler et al.l 
(|2001bD . the A6437 A feature is significantly blended with 
a silicon feature at A6437.71 A. We have attempted to 
reduce the impact of this blend and recover reliable Eu 



abundances from this feature in the following way. We 
measured the Eu abundance for Arcturus from the Eu II 
A6645 line and set that as the "true" Eu abundance: log 
e(Eu) = 0.25. We then generated a synthetic spectrum of 
the A6437 feature using Arcturus's atmospheric param- 
eters, the Eu abundance from the A6645 line, and the 
Si abundance we have adopted for Arcturus (and which 
served as the basis for our astrophysical log gf's; we also 
used an E.P. of 5.863 eV, retrieved from VALlQ. We 
then altered the log gf of the A6437.71 Si feature until the 
synthetic spectrum matched the Arcturus spectrum, and 
used this log gf (—2.30) for the Si feature in the synthesis 
line list for Eu A6437. Lastly, in the spectrum synthe- 
sis analysis of our sample stars, we set the Si abundance 
of this blend feature to match that for each star as de- 
termined from the EW analysis of other Si lines in its 
spectrum. 

As for La, we verified our method of analysis with 
the S un. Spectrum synt hesis using the solar flux atlas 
from iHinkle et all (|2000D resulted in log e(Eu) = 0.53 
and 0.49 for A6437 and A6645, res pectively. These val - 
ues are consistent to those found bv lLawler et al.l (|2001bh 
for the same lines: 0.55 and 0.54. For Arcturus, our 
[E u/Fe] = +0.23 i s in g ood a greement wi t h tha t found 
by iJohnson et al.l (|2012f ) and lYong et al.l (|2005f) . again 
using similar analysis techniques ([Eu/Fc] = +0.29 in 
both studies, from only the A6645 line). Our value of 
[Eu/ H] = —0.27 f o r Arc turus is also in good agreement 
with iSmith et al.l (pool , who found [Eu/H] = -0.30. 
Table [7] presents the [Eu/Fe] ratios for our sample as 
determined from spectrum synthesis. Again, the abun- 
dance uncertainties associated with each line represents 
the measurement uncertainty from the synthesis tech- 
nique. We note that the S/N ratios of the spectra of all 
four NGC 1245 stars were too low to determine reliable 
abundances for the relatively weak (<20 mA in these 
stars) Eu II features. Mean cluster [Eu/Fe] abundances, 

12 See http ://vald.astro.univie.ac.at/~vald/php/vald.php 
IKupka et all l|20UUi ). 
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Table 9 

Uncertainties in abundance ratios due to atmospheric parameters 







T+100 K 


log g+0.2 dex 


[M/HJ+0.2 dex 
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El 




Aabund. 


Aabund. 


Aabund. 


Aabund. 


N6939 31 


Fe/H] 
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01 
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with (1 a) standard deviations, are shown in Table [SJ 
(As for La, the 1 a value for clusters based on the anal- 
ysis of a single star represents the uncertainty (standard 
deviation) in its abundance from Table [7]) The bottom 
panel of Figure [5] shows an example synthesis of the Eu 
II 6645 feature. 

It can be seen in Table [7] that the Eu abundances de- 
termined from the A6437 feature generally agree well (i.e, 
with differences comparable to or smaller than the mea- 
surement uncertainties) with that of the A6645 feature 
for most stars. However, larger discrepancies are evi- 
dent for some stars, especially those in clusters NGC 188 
and NGC 7142. While differing line abundances could 
be due to uncertainties in continuum placement or less- 
than-optimal S/N ratios for some stars, the stars exhibit- 
ing the largest abundance discrepancies are generally the 
more metal rich stars in the sample, with [Fe/H] values 



of ^0.10 or higher. For such metal-rich stars, our fix for 
neutralizing the Si blend in the A6437 feature may break 
down. 

Eu abundances were determined using measured EWs 
and the blends driver in MOOG, similar to lanthanum. 
As before, differences between EW and synthesis abun- 
dances were plotted against stellar [Fe/H], T e gf and log 
g. Both Eu features show a trend of increasing difference 
with increasing [Fe/H], but no trend with T e ff or log g 
(maybe a slight trend with log g for A6437). The mean 
difference for both lines is ~0.1 dex, with the smallest 
abundance difference for metal poor stars (^0 dex) and 
the largest abundance difference (~0.2 dex) for the more 
metal-rich. 



3.3. Error analysis 



12 Jacobson & Friel. 



Table 10 

Cluster abundances from the literature 



Cluster 


Ref 


[Zr/Fe] 


[Ba/Fe] 


[La/Fe] 


[Eu/Fe] 


Notes 


Be 18 


Yong et al. (2012) 


—0 21 +0 1 3 


+0 30+0 06 


+0 34+0 08 


+0 30+0 06 


2 stars 


Be 21 


Yong et al. (2012) 


— 0.11±0.08 


+0.59+0.01 


+0.56+0.01 


+0.31+0.08 


2 stars 


Be 22 


Yong et al. (2012) 


-0.14 


+0.61+0.05 


+0.38+0.02 


+0.26+0.10 


2 stars 


Be 32 


Yong et al. (2012) 


-0.06±0.07 


+0.29±0.12 


+0.44+0.02 


+0.31+0.11 


2 stars 


Be 32 


D'Orazi et al. (2009) 




+0.24±0.15 






9 stars 


Be 32 


Carrera & Pancino (2011) 




+0.51±0.12 


-0.14+0.07 




2 stars 


M 67 


Yong et al. (2005) 


-0.28±0.03 


-0.02±0.04 


+0.11+0.02 


+0.06+0.03 


3 stars 


M 67 


Tautvaisiene et al. (2000) 


-0.17+0.09 


+0.07+0.09 


+0.13+0.09 


+0.07+0.09 


9 stars 


M 67 


D'Orazi et al. (2009) 




+0.04±0.05 






10 stars 


M 67 


Maiorca et al. (2011) 


+0.04+0.05 




+0.06+0.04 




10 stars 


M 67 


Pancino et al. (2010) 




+0.25±0.06 


+0.05+0.06 




3 stars 


NGC 1817 


Reddy et al. (2012) 


+0.08±0.05 


+0.13 


+0.12+0.03 


+0.13 


3 stars 


NGC 2141 


Yong et al. (2005) 


+0.63+0.06 


+0.91 


+0.57+0.07 


+0.17 


1 star 


PWM 4 


Yong et al. (2012) 


+0.11+0.06 


+0.46 


+0.22+0.07 


+0.12 


1 star 



As in our previous work, we estimate the uncertainties 
in atmospheric parameters to be 100 K in T e ff, 0.2 dex 
in log g, and 0.2 km s _1 in microturbulent velocity. We 
have calculated uncertainties in [Fe/H] and [X/Fe] ratios 
due to these uncertainties in atmospheric parameters for 
six stars in our sample that span a wide range in at- 
mospheric parameters and metallicities. The results are 
presented in Tabled We did not consider the covariance 
terms in the error analysis, rather we treated the uncer- 
tainty an each parameter as though they were indepen- 
dent of one another. The uncertainty calculations for La 
and Eu lines were performed using EWs, not spectrum 
synthesis. We include here also the abundance changes 
due to an uncertainty of 0.2 dex in the model atmosphere 
[M/H] value. This is meant to be indicative only of the 
impact that such an uncertainty would have; given that 
we always matched the [M/H] of the model to the iron 
abundance of each star, this uncertainty is overly conser- 
vative. As can be seen, [X/Fe] uncertainties are generally 
smaller than 0.10 dex, but occasionally as large as 0.20 
dex. 

To estimate the uncertainty in the Eu abundance cal- 
culated from the A6437 line due to the Si blend, we gen- 
erated a synthetic spectrum for each of the stars in Table 
[9] with the Si abundance varied by the line-by-line dis- 
persion in its Si abundance as measured from EWs of Si 
lines. The best fit Eu abundance was then compared to 
the abundances given in Table[7l Line-by-line dispersions 
in Si abundances ranged from 0.08 to 0.16 dex. Corre- 
sponding changes in the Eu abundance ranged from to 
^0.15 dex, with a median change of only 0.02 dex. 

3.4. Anomalous stars 

A few clusters contain stars that exhibit strikingly dif- 
ferent n-capture abundance patterns. 

Be 22: Star 643 is much more enhanced in [Zr/Fe] and 
[Ba/Fe] than star 414. In fact, the abundances of several 
other elements in 643 differ significantly from those of 
414, and the general line-by-line dispersion in element 
abu ndanc es is also higher (even for iron; we also note 
that lY12l found similar behavior). Given that 643 is high 
up on t he red giant branch relative to 414 (e.g., Figures 
we have taken the abundances of star 414 
to be more reliable, and therefore we have adopted the 
abundances of star 414 as representative of those for Be 
22 as a whole. The abundances reported for Be 22 in 



Table [8] are exactly those of star 414. 

NGC 2141: Star 1348 shows systematically larger en- 
hancements in [Zr/Fe], [Ba/Fe] and [La/Fe] than star 
1007, though their [Eu/Fe] ratios are nicely consistent 
(see Figure [1} . They have very similar atmospheric pa- 
rameters, met allicity, and other element abundances. We 
also note that I Yong et all (|2005l ) also found high [Zr/Fe], 
[Ba/Fe], and [La/Fe] ratios for star 1348. Given the sim- 
ilar evolutionary states of these two stars, we do not 
choose one and exclude the other as we have done for 
Be 22. Instead, we have averaged the abundances of the 
two stars together, and therefore the mean abundances 
of this cluster show large dispersions in the subsequent 
discussion. 

NGC 6939: Stars in this cluster show a pattern in 
which one star, 190, is consistently enhanced in [X/Fe] 
ratios, and another (121) is consistently under- abundant 
relative to the remaining two stars. Star 190 resides in 
the clump and shows [s/Fe] ratios at least 0.4 dex larger 
than the other stars. As a result we do not include it in 
cluster mean abundance calculations. Star 121 is more 
of a puzzle, as it is not systematically different from the 
other two cluster members for the other elements. How- 
ever, we note that the spectroscopic T ff that we found 
by minimizing Fe I abundances versus line E.P. is 300 
K cooler than the photometric T c g (Jacobson et al., in 
prep.), while the photometric and spectroscopic temper- 
atures for the other stars agree within 50 K. Therefore, 
it is possible that star 121's discrepent [X/Fe] ratios are 
due to offsets in its atmospheric parameters. 

Other stars show less pathological differences from 
their counterparts. For example, star 398 in NGC 2355 
is 0.3 dex more enhanced in [Ba/Fe], while star 1224 in 
NGC 188 also shows larger (~0.2-0.3 dex) enhancements 
in [Ba/Fe] and [La/Fe] than the other NGC 188 stars, 
even after taking into account line-to-line abundance dis- 
persions. The four stars in NGC 1245 show large varia- 
tions in [La/Fe], but this may be the result of relatively 
low S/N ratios for the spectra (indeed, they were deemed 
too poor to reliably measure Eu abundances). In general, 
star-to-star variations in abundances of the r-process el- 
ement Eu are much smaller than for the s-process ele- 
ments. For example, the pairs of stars in Be 22 and NGC 
2141 that differ in [s/Fe] ratios have similar [Eu/Fe] ra- 
tios, and the Eu abundance for star 190 in NGC 6939 
is consistent within the uncertainties of the other stars 
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Figure 3. [Ba/Fe] (top panel), [La/Fe] (middle panel), and [Zr/Fe] (bottom panel) ratios for our cluster sample (filled squares) as a 
function of cluster age. The open blue circles represent clusters Be 20 and Be 29 from Yong et al. (2005) shifted on to our abundance scale 
(see text for details). The equations for least squares fits to the data are indicated in all panels, along with the R- values, p- values, and 
errors of the fit. 

in the cluster. Instead, the two stars in Be 21 differ in 
[Eu/Fe] by 0.2 dex, but at a 2-er level. Not surprisingly, 
the star-to-star scatter is larger for NGC 188 and NGC 
2355 stars, with standard deviations in [Eu/Fe] of 0.12 
and 0.28 dex, respectively. 



3.5. Comparison to other studies 

To our knowledge, abundances of neutron capture el- 
ements are available in the literature for eight clusters 
in our study. Table [TO] gives a summary of these liter- 
ature results. Not surprisingly, M 67 is the most well- 
studied with two sets of measures for Zr, Ba , La, and Eu 
(jYong et all 120051: iTautvaisiene et all |2000j), and three 
additional studies o f Ba and La abundances (D09; Mil; 
iPancino et alJl2010( ). For this cluster, our [Ba/Fe] ratio 
is consistent with the roughly scaled-solar values found 
by these studi es, save for the [Ba /Fe] = +0.25±0.06 
value found by IPancino et al.l ((2010). Our [La/Fe] ratio 
of —0.15 ± 0.04 for M 67 is systematically lower than 
those of other studies, which found [La/Fe] ^0.05 to 
0.13 dex. Our [Eu/Fe] measures for M 67 stars are also 
lower than the find ings of other studies in t he literature 
( e.g., +0.07± 0.09- lT^vaisiene et al.ll2000l . +0.06±0.03 
- I Yong et alll2005l) . As for Zr, our value is consistent 
with that of IM 111 but is larger than those values found 
by Yong et al. and Tautvaisiene et al. 

The next best-studied cluster is Be 32, with lit- 



erature abundance s agai n from El IPOflL [Mill , and 
iCarrera fc Pancinol (|2011h . As for M 67, our mean clus- 
ter [Ba/Fe] ratio agrees well with those of IY12I and ID09I 
(within 0.05 dex), but is ~0.25 dex lower than th at of 
Carrera & Pancino. Only Carrera & Pancino and IY12I 
presented [La/Fe] ratios for Be 32; our value is in good 
agreement with the former (within 0.06 dex), but is lower 
than the latter by 0.5 dex! S imila rly our [Eu/Fe] ratio 
is 0.3 dex lower than that of IY12L but our scaled-solar 
[Zr/Fe] is consistent with their findings. 

The recent work by iReddv et al.l (|2012D also presented 
neutron capture abundances for NGC 1817. Here we find 
consistent [La/Fe] ratios (within 0.02 dex) for the cluster, 
while our [Ba/Fe] and [Zr/Fe] ratios are consistent within 
the uncertainties. Our Eu abundances, however, are very 
different: their [Eu/Fe] = +0.13 is 0.26 dex larger than 
our value. 

In addition to the clusters mentioned above, clusters 
NGC 2141, Be 18, Be 21, Be 22 and PWM 4 are in 
common with Yong et al. (2005, 2012) - recall the stars 
and spectra for the three Berkeley clusters and star 1348 
in NGC 2141 are those of Yong et al., and we have in 
addition two M 67 stars in common (105 and 141). A 
star by star comparison of the 12 stars common be- 
tween our two samples shows that differences between 
our abundances (in the sense This Study — Yong) for 
[Ba/Fe], [La/Fe], [Eu/Fe] and [Zr/Fe] are +0.04±0.19, 



14 



Jacobson & Friel. 



1 



0.5 



a. 

\ 

CO 
CQ 



-0.5 



D'Orazi dwarfs 
D'Orazi giants 
Yong 

This Study 
Lit. 




10 8 



10» 

log(age) (Gyr) 



10'° 



Figure 4. [Ba/Fe] versus log(age) for our clusters (open squares), those of Yong et al. (2005, 2012; red triangles), and the sample of ID09I 
(dwarfs - filled blue circles, giants - open blue circles). Other clusters from the literature are represented by small crosses. See text for 
more information. 



-0.29±0.12, -0.16±0.13, and +0.26±0.10 dex, respec- 
tively. The standard deviations of this difference for La, 
Eu, and Zr are comparable to individual star uncertain- 
ties, while the scatter for Ba is larger, indicating that the 
differences between our studies are systematic for La, Eu 
and Zr, but not for Ba. For the discussions that follow, 
we have made use of these systematic differences to place 
two additional clusters studied by Yong et al. (2005), Be 
20 and Be 29, on to our cluster abundance scale. 

4. DISCUSSION 

In order to minimize possible systematic uncertainties 
in the cluster abundance results as much as possible, we 
limit our discussion to the results of this wor k with the 
addition of Be 20 and Be 29 from lYong et all (|2005f ). as 
just discussed. For the first discussion, that of abundance 
trends with age and R oc , our sample size is comparable 
to those of ID09I and IM11I and yet is almost completely 
distinct from them, with only two clusters in common 
to each (M 67 and Be 32). Therefore, this dataset pro- 
vides an excellent opportunity to look for the anticorre- 
lation of cluster s-pro cess abun dance ([s/Fe]) with clus- 
ter age found b y ID09I and iMlll . and further explored in, 
iMaiorca et~aT] (|2012fl . 

4.1. Cluster n-capture abundance trends with age 

The first s-process element shown to have enhance- 
ments anti-correlated with open cluster age was barium 
(D09), so we begin with it here. The top panel of Figure 
[3] shows mean cluster [Ba/Fe] ratio s from Table EQ plot - 
ted against cluster age (taken from lSalaris et ah 12004). 
In this figure, our sample is represented with squares, 



while open blue circles represent the outer disk clusters 
Be 20 and Be 29 from Yong et al. (2005), shifted to our 
abundance scale. A simple linear regression analysis on 
the distribution of the points resulted in the equations 
for lines of best fit through the data also printed in each 
panel. Also included are the R-value, standard error of 
the estimate for regression (err; i.e., the degree of scatter 
of the points around the regression line), and the two- 
tailed probability estimate (p- value), a measure of the 
statistical significance of the slope. Conventionally, p- 
values smaller than 0.05 or 0.01 are taken to be rejections 
of the null hypothesis, which in this case is that there is 
no trend in [X/Fe] ratio with cluster age. As can be seen 
in Figure [3l our data do show a statistically significant 
trend of increasing [Ba/Fe] wit h decreasing cluster age, 
as shown in the sample of ID09L If the oldest cluster, Be 
17, is excluded, the p-value increases from p=0.010 to p 
= 0.048, decreasing the significance of the trend. 

Re sults from our sample also agree with those of 
IY12L who found a similar slope of [Ba/Fe] with age, of 
—0.03 ± 0.01, from their sample of 10 clusters combined 
with another 22 measurements taken from the literature. 
As they note, this trend is not as pronounced as that 
found in D'Orazi et al. The advantage of the work pre- 
sented here is that it includes 3 more clusters with ages 
greater than 6 Gyr. Notably, our [Ba/Fe] ratios for these 
old clusters are all roughly solar, reinforcing the signifi- 
cance of the trend of decreasing Ba with cluster age. Nev- 
ertheless, the correlation shows substantial scatter about 
the mean trend. This can a lso be seen in Figure|3J which 
replicates Figure 1 in lD09l (open and filled blue circles), 
with our sample (open squares) also included. Here, we 
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Figure 5. [La/Fe] (top) and [Zr/Fe] (bottom) abundance ratios versus age for our sample (open squares), those of Yong et al. (2005, 2012; 
red triangles) shifted to our abundance scale , and those of Mil (blue open circles). Error bars have been omitted for clarity. As described 
in the text, the [s/Fe] for clusters from ,\l I lj fit nicely within the distribution of our sample. 



include the full cluster sample of Yong et al. (2005, 2012; 
red triangles ), alo ng with clusters from the literature as 
compiled in IY12I ( their Table 13) represented by small 
crosses. We follow IY121 here in that we treat each mea- 
sure of an individual cluster independently, and therefore 
some clusters are represented by multiple points in this 
figure. It can be seen that although each sample spans 
a slightly different range in age, they follow the same 
general trend of increasing [Ba/Fe] with decreasing age. 
This is even seen in the inhomogeneous cluster sample 
drawn from the literature. 

The middle panel of Figure [3] shows cluster [La/Fe] ver- 
sus age for our sample, again with two outer disk clus- 
ters fr om Y ong et al. (2005) adjusted to our abundance 
scale. IM11I found a trend of increasing [La/Fe] with de- 
creasing cluster age, but this was based on only six clus- 
ters in their sample for which they could determine La 
abundances (the trend of [s/Fe] with cluster age is much 
more convincing for the s-process elements Ce and Y, for 
which they have abundances for all 19 clusters in their 
sample). We note that the clusters for which they de- 
termined La abundances covered a much more limited 
age range (and a primarily younger one) than that pre- 
sented here. The larger sample of Figure |3] clearly shows 
no trend of [La/Fe] with cluster age: the p-value is quite 
large. Indeed, while Be 17 shows a subsolar [La/Fe] ratio 
that drives some slope to the data, the rest of the clusters 



show an essentially flat distribution, again with a scatter 
o f roug hly 0.3 dex. 

IY12L on the contrary, found a positive correlation of 
[La/Fe] with age, with a slope of 0.05 ± 0.01, based on 
their determinations for 10 clusters and literature values 
for 6 more. However, r ecall that our values of [La/Fe] for 
stars in common to the lY12l analysis showed a syste matic 
offset of 0.29 dex, ours being lower. When the lY12l values 
are adjusted to our abundance scale, the trend with age 
disappears entirely, yielding results consistent with ours, 
as can be seen in the top panel of Figure [3 Here, the 
sample of IY12I . shifted to our abundance scale , is in di- 
cated by red triangles, while the six clusters in IM11I are 
blue circles We conclude that none of the av ailabl e data 
from these larger samples (that is, ours and IY12[ ) show 
evidence for age trends in La abundance, once systematic 
differences are taken into account, and the distribution 
with age of the IM11I sample is consistent with this obser- 
vation. 

IM11I also claimed a trend in their cluster [Zr/Fe] ratios 
versus cluster age, but as for La this was based on a sub- 
set of their sample of only 6 primarily younger clusters. 
We see no trend in our cluster [Zr/Fe] ratios versus age; 
indeed, the distribution is very similar to that of [La/Fe] 
(bottom panel of Figure |3]) . Interestingly, the dispersion 
of [Zr/Fe] ratios with age is fairly tight at 0.1-0.2 dex for 
clusters around 6-7 Gyr of age, with the dispersion in- 
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Figure 6. Top panel: [Eu/Fe] versus age; symbols same as in Figure [3] Bottom panel: [Eu/Fe] for our cluster sample (filled squares), 
along with those of Yong et al. (2005, 2012; red triangles) shifted to our abundance scale, and other clusters from the literature (crosses). 



creasing with decreasing cluster age. We note the cluster 
with [Ba/Fe] -0.70 with the huge errorbars is NGC 2141, 
whose large and discrepant Zr abundances have already 
b een d iscussed. 

IY12I also found a slight anticorrelation of [Zr/Fe] with 
cluster age, with a slope of —0.03 ± 0.01, based again on 
their sample plus results from the literature, the majority 
of which were taken from our previous papers. We note 
that these earlier Zr abundances have been re vised and 
updated in this paper, as discussed in Section 3. lY12l also 
noted that the tr end of Zr with age was more moderate 
than that cited by I Mill H owever, as with La, the [Zr/Fe] 
abundances for the iMllI sample are consistent with the 
distribution of points shown in the bottom panel of Fig- 
ure [5J and support the conclusion that there is no trend 
of [Zr/Fe] abundances with age among the open clusters 
from this larger sample. 

Very few determinations of Eu exist for stars in open 
clusters, yet the behavior of this r-process element rela- 
tive to the s-process elements we have considered thus far 
is especially interesting. Figure |6] shows cluster [Eu/Fe] 
trends versus age for our sample. Our data show a trend 
of increasing [Eu/Fe] with increasing cluster age but only 
at a marginally statistically significant level. What is 
more notable in these figures is how the dispersion in 
cluster [Eu/Fe] ratio varies with age: for example, the 
dispersion in [Eu/Fe] among 4 Gyr old clusters spans 
nearly 0.6 dex, while the dispersion of the older and 



younger clusters is 0.3-0.4 dex. Clearly, this could just 
be a result of the small size of our cluster sample, but 
it would be interesting to see if thes e disp ersions hold in 
larger homogeneous cluster studies. IY12I found no trend 
of [Eu/Fe] with age in their analysis, based on their re- 
sults for 10 clusters plus four others from the literature. 
Those results are consistent with our larger sample, as 
can be seen in the bottom panel of Figure [6l 

As mentioned in the introduction, Ida Silva et al.1 
(|2012T) found trends of increasing [s/Fe] with decreas- 
ing age for dwarf stars in the solar neighborhood, most 
notably for Ba, but also for Zr, Y and Ce. Of these el- 
ements, only Y shows this anticorrelation across the full 
age range of their sample; the anticorrelation is present 
for the other elements only for the stars you nger t han 
the Sun (4.5 Gyr). They argue, as pi and ImIH do, 
that this indicates that the s-process yields of AGB stars 
has changed over time (as well as [Fe/H] of the AGB 
stars), and the fact that Ba shows the strongest trend 
with age implies that this change in yields affects the 
heavy s-process elements more than it does the light 
ones. While we find the c onsis tency betw een the re- 
sults of Ida Silva et all (|2012D and iDOiMMTTl compelling, 
we would argue that the [La/Fe] trends seen in Figures 
[3] and [5] of this work may indicate that the source of an 
[s/Fe] abundance trend with age does not affect all heavy 
s-process elements equally. 

Indeed, the picture becomes more complicated when 
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Figure 7. [Ba/Eu] (top panel). [La/Eu] (middle panel), and [Zr/Eu] (bottom panel) versus age. Symbols are the same as in Figure|3] 



other (larger) solar neighborhood dwarf star studies are 
considered. We performed a regression analysis of [X/Fe] 
versus age for sola r neighborhood (thi n dis k) dwarf stars 
from t he studies of lReddv et al.l (|2003D and lBensbv et al.l 
( 2005), considering the full age range of the ir samples as 
well a s stars 4.5 Gyr or younger. For the IReddv et al.l 
(2003) sample, [Zr/Fe] showed a statistically significant 
trend with decreasing stellar age for both age ranges con- 
sidered, while [Nd/Fe 13 showed a statistically significant 
trend with (increasing) age only wh en the full age range 
the sample was considered. For the lBensbv et al.l (2005) 
sample, only [Ba/Fe] showed a statistically significant 
trend with (decreasing) age, for both age ranges con- 
sidered. None of the ot her s-process specie s studied by 
th ese two group s (Y for iBmsbv et al.l [2005 Ba and Ce 
for IReddv" et al. 2003|) showed trends with age. (The lack 
of a trend of [Ba/F e] with decreasing cluster age for the 
IReddv eta l. (2003) sample may be due to the fact that 
that sample contained very few stars with young (<1- 
2 Gyr) ages. If the distribution of [Ba/Fe] ratios with 
age is non-linear, and in fact only becomes steeply anti- 
correlated with age for t he youngest objects (e.g., the 
youngest clusters in ID09D , then such a trend would not 
appear in a sample of older stars.) 

Regarding r -proc ess abundance trends with age, 
Ida Silva et al.l ()2012l ) did not determine Eu abundances 

13 Technically, Nd is ~50% s-process and ~50% r-process 
dBurris et al.ll20utjf) . 



for their sample. However, they did report abun- 
dances for Sm, which is also representative of the r- 
process. They also reported a trend of increasing 
[Sm/Fe] with increasing age for their sample, consistent 
with the Eu result s for our open cluster sample. The 
IReddv et al.l (|2003D sample also showed statistically sig- 
nifi cant trends of [ Eu/F e] with increasing stellar age, but 
the iBensbv et al.l ([2005) sample did not. 

In Figure [7] we explore trends of [X/Eu] ratios versus 
age for Ba, La and Zr individually. Trends of increas- 
ing [Ba/Eu] and [La/Eu] with decreasing cluster age can 
be seen, while that for [Zr/Eu] is less convincing. These 
trends are con s istent with the [s / Sm] versus age trends in 
Ida Silva et al.l (|2012l ) , although the trends in their sample 
are especially pronounced for stars younger than the Sun. 
We interpret this as evidence of the increasing impor- 
tance/efficiency of s-process enrichment relative to the 
r-process in the Milky Way thin disk over time. 

This analysis did not take into account uncertainties 
in the ages of the clusters, which can sometimes be sig- 
nificant. W e have chosen to us e the cluster ages as de- 
termined in lSalaris et al.l (|2004D which r est on a calibra- 
tion o f the morphological age indicator of Janes & Phelps 
(1994). While this method may not be as accurate for 
any given cluster as direct fitting of isochrones to clus- 
ter color magnitude diagrams, it has the advantage of 
placing all the cluster ages on a uniform relati ve scale. 
The typical age uncertainty of the ISalaris et al.l calibra- 
tion is 15%, which for our clusters is less than 1 Gyr. 
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To give an indication of the impact of age uncertainties 
on the results, we varied the age of the oldest cluster, 
Be 17, by it s uncertainty of 2.77 Gyr as given by the 
ISalaris et alj calilbration. This estimate of the uncer- 
tainty in age offers an extreme case, both in terms of the 
ISalaris et al. | calibration a nd direct age determinations 
(e.g-.TBragagli a et al.ll2006l found an age of 8.5 - 9 Gyr for 
Be 17 from isochrone fits). Making Be 17's age younger 
or older both served to slightly decrease the statistical 
significance of the trends shown in Figure 3, e.g., p=0.013 
and 0.017 for [Ba/Fe]. Coupled with the leverage Be 17's 
large age gives to the fitting of element abundance trends 
with age, this example provides an indication of the ro- 
bustness of such trends to age erro rs. Adopting t he more 
typical 15% age uncertainty of the ISalaris et al.l calibra- 
tion would have significantly less impact on any of these 
correlations. 

Lastly, we also performed a regression analysis of our 
sample's [X/Felratios as a function of log age, such as 
seen in Figure 21 to assess how sensitive the statistical 
significance of the trends were to how age was repre- 
sented. We found that, in general, the p-values increased 
for each element considered here, decreasing the statisti- 
cal significance of each trend. For example, the p-value 
of [Ba/Fe] versus log age was 0.041, as opposed to 0.010 
as seen in Figure [3] However, the opposite is found in 
an analysis of the ID09I sample: the p-value of their clus- 



ter [Ba/Fe] ratios versus log age was 0.0003, compared 
to p = 0.011 found versus age. This i s no doubt due to 
the much larger age range of the ID09I sample relative to 
ours (the youngest being 35 Myr old, as opposed to our 
700 Myr). These changes in p-value do not change the 
primarily conclusions described above. 

4.2. Cluster n-capture abundance trends with R gc 

The open clusters have often been used as tracers of 
the elemental abundance distributions in the Milky Way 
disk. In terms of the traditional characteristics used to 
categorize disk populations, the properties of the open 
clusters show overall consistency with those of the thin 
disk (Friel 1995). Their kinematics, and in particular the 
systemic rotational properties of the open cluster popu- 
lation, are very similar to those of the thin disk field star s 
and other disk tracers (jHronl Il987t iScott et all 119951 ). 
Cluster ages, spatial distribution, and concentration to 
the Galactic plane a re thin disk-like. As demonstrated by 
many papers (e.g., Uacobson et al J 120111 : iMagrini et al.l 
120091 the elemental abundance ratios of open clusters 
follow the trends shown by field stars of the thin disk. 
Although there may be a few open clusters whose in- 
dividual properties dist inguish themselves in interesting 
ways (e.g. NGC 6791 - IPlatais et aTllMlt iGeisler et al.1 
l2012f ). the evidence for the open cluster population as a 
whole is that it provides an important and useful sam- 
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pie with which to study the evolution of the thin disk 
properties. 

It is established by now that open clusters (and other 
disk populations) show that the thin disk's mctallic- 
ity gradient, which decreases roughly with a slope of 
—0.07 dex kpc" 1 throug h the solar neighborhood (e.g., 
iCarrera fc Pancinol |20Tlj ) . flattens out at [Fe/H] —0.3 
around R gc ~13 kpc. Cluster [X/Fe] ratios are essentially 
flat across the full R gc range spanned by open clusters for 
the light, a, and Fe-peak elements, indicating that they 
all follow Fe fairly closely. It is therefore interesting, es- 
pecially in light of the trends of at least some of them 
with cluster age, to investigate the trends of n-cap ture 
element abundances with R gc . To this end, IM11I and 
Yong et al. (2005, 201 2) ha ve all investigated n-capture 
abundances with R gc . IM11I remarked that the s-process 
elements they studied have distributions similar to the 
other element groups studied (i.e., a, Fe-peak) - that 
is, that cluster [s/Fe] is essentially independent of R gc . 
That said, they note that the trends with R gc for young 
and old clusters are different, as expected based on the 
abundance trends with age found in their sample (see, 
e. g., th eir Figures 6 and 7 for Ce and Y abundances). 

Y12I investigated radial abundance trends in clusters 



with a focus on distinguishing behavior in the solar neigh- 
borhood from that in the outer disk (taken to be Rgc > 
13 kpc). For the neutron-capture elements they find that 
[Zr/Fe] shows no trend with R gc , but for Ba, La, and Eu, 
there is a suggestion that the outer disk clusters follow a 



different gradient from that of the local sample of clus- 
ters. In particular, for La and Eu, the outer disk clusters 
appear enhanced relative to those inside ~ 10 - 13 kpc. 

We present [Ba/Fe], [La/Fe] and [Zr/Fe] ratios as a 
function of cluster R gc for our sample in Figure [5] In 
each panel, we distinguish older (4 Gyr and greater) and 
younger (younger than 4 Gyr) clusters by open and filled 
symbols, respectively, though the details of regression 
analyses given in each panel is for the entire cluster sam- 
ple. Given the narrow R gc range of the bulk of our clus- 
ter sample (R gc ^8-14 kpc), we do not draw conclusions 
about general abundance distributions in the disk as a 
whole. Rather, we present a linear regression analysis 
of our original cluster sample alone (with R gc <14 kpc) 
and include the results in each panel for comparison to 
the gradients presented by IY12I for what they term the 
"inner disk" (R gc <13 kpc; their Figure 24). Generally, 
the gradients we fin d for [La/Fe] and [Zr/Fe] are larger 
than those found bv lY12l while our slope for [Ba/Fe] is 
shallower than theirs. None of these slopes are statisti- 
cally significant. Figure [9] sh ows [L a/ Fe] an d [Zr/Fe] for 
our sample along with that of lY12l and IM 111 as a function 
of R gc . Again, the distributions of all three samples are 
consistent with one another. 

Figure ITD1 shows [Eu/Fe] versus R gc , again with the 
clusters divided into you nger (filled symbols) and older 
(open symbols) age bins. IY12I reported a possible trend 
of enhanced [Eu/Fe] in the outer disk. However, with 
their [Eu/Fe] ratios adjusted to our abundance scale, the 
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four clusters beyond R gc ~13 kpc in their sample have 
abundance ratios consistent with the inner disk clusters 
(bottom panel of Figure [TOl) . Therefore, in comparison 
to the conclusions from IY12I . we note that the system- 
atic differences between our abundance scales for [La/Fe] 
(-0.29 dex) and [Eu/Fe] (-0.16 dex), when taken into 
account, would serve largely to remove the differences 
they see between the inner and the outer disk abun- 
dances, yielding trends consistent with those seen in Fig- 
ures M and [H 

That said, the younger and older populations of clus- 
ters appear distinct in [Eu/Fe] versus R gc (Figure fTUj) . 
The younger clusters (filled symbols) do not exhibit sta- 
tistically significant trend with R gc , but the trend with 
R gc for the older clusters (open symbols) is marginally 
significant (p=0.017). We found, however, that these 
trends are largely driven by the radial variation of clus- 
ter [Fe/H] rather than of [Eu/H], indicative of a possible 
flattening of the thin disk abundance gradient over time. 
It would be interesting to explore possible variations in 
the radial distributions of younger and older clusters in 
larger, homogeneous cluster samples that better probe 
the full Rg C and age range of the Galactic disk. 

4.3. Other abundance ratio trends 

As discussed in earlier sections, barium and lanthanum 
are both heavy s-process elements, so it seems reasonable 



to expect their abundances to scale together. Figure [TT] 
shows [Ba/Fe] versus [La/Fe] (top panels) for our clus- 
ter sample. Here again, we distinguish younger clusters 
(<4 Gyr) from older clusters (>4 Gyr) with filled and 
open symbols, respectively. It can be seen from the top 
left panel that [Ba/Fe] and [La/Fe] do appear to scale 
together, though with a constant offset from a 1-to-l 
correlation. [Zr/Fe] also seems to be well correlated with 
[La/Fe] for most clusters in our sample as well (top right 
panel). Also, younger clusters in our sample are slightly 
shifted to higher [Ba/Fe] ratios from the older clusters as 
a result of the [Ba/Fe] correlation with age we have al- 
ready discussed, but with much overlap (especially when 
errorbars are taken into account). 

As Zr is a light s-process element, the ratios of [Ba/Zr] 
and [La/Zr] provide information about the efficiency of 
s-process element production in the AGB stars that en- 
riched the gas from which clusters in the disk formecQ 
The bottom panels of Figure [11] shows the ratios of Ba 
and La to Zr as a function of cluster [Fc/H] for our sam- 
ple with the clusters distinguished by age. Neither age 

14 As mentioned previously, as a light s-process element, Zr can 
also be produced via the weak s-process in mass ive stars (e.g., 
Busso & Gallino 1985). However, as mentioned in Major ca et al.l 
120121) . the estimated weak s-process contribution to the solar sys- 
tem l ight element abundance distribution is small (Scrminato et al. 
[20091) . 
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group shows a statistically significant trend of [X/Zr] 
with [Fe/H]. Indeed, for the cluster sample as a whole, 
these [hs/ls] versus [Fe/H] trends are rathe r flat, similar 
to those found for thin disk dwarf stars in iReddv et al.l 
(2003; their Figur e 14) We note that plots of [Ce/Y] ver- 
sus [Fe/H] for the lMll j sample also show a flat distribu- 
tion with metallicitv. IReddv et al.l say, "The similarity 
of the slopes of [X/Fe] versus [Fe/H] for light (Sr, Y, Zr) 
and heavy (Ba, Ce) elements may reflect either unchang- 
ing relative contributions from AGB stars or a change in 
these contributions that is offset by a change in the weak 
s-process contributions from the massive stars. On the 
assumption that the AGB stars are the controlling influ- 
ence, the unchanging abundance ratio of heavy to light 
elements indicates that the exposures to neutrons in the 
s-process site is essentially independent of the metallicity 
of the contributing AGB stars." If this interpretation is 
true, it implies that s-process yields of AGB stars do not 
change as a function of [Fe/H] in the metallicity range 
spanned by our sample, but may still vary as a function of 
AGB star age (i.e., mass), as mentioned in Ida Silva et al.1 
(|2012f) . This interpretation i s con sistent with the model 
proposed by iMaiorca et al.l (|2012f ). but again the situa- 
tion is confused by the uncertainties in the abundance 
trends with age seen for the heavy and light s-process 
elements among the various studies. 

Figure [H plots [Ba/Fe], [La/Fe] and [Zr/Fe] versus 
[Eu/Fe] for clusters in this study. In addition, we also 



show [La/Fe] versus [E u/Fe] for red horiz ontal branch 
stars from the study of lAfsar et al.l (|2012[ ). Plots of s- 
process elements versus r-process elements reveals the 
relative contributions of the s- and r-processes to the 
chemical enrichment of the Milky Way disk. Another 
presentation of [s/r] ratios relative to [Fe/H] is given in 
Figure HH As can be seen in the bottom panel, the older 
clusters in our sample (open squares) have scaled solar 
[La/Eu] ratios, except for the two clusters with super- 
solar metallicity (NGC 7142 and NGC 188). The younger 
clusters, on the other hand, generally show enhanced 
[La/Eu] ratios. A similar pattern is shown for [Ba/Eu], 
too, shifted relative to that for La (recall the top left 
panel of Figure [TT|) . The distribution of [La/Eu] ver- 
sus [Fe/H] i s similar to that shown by the RHB stars of 
lAfsar et al.l as well, indicating that metal-rich RHB stars 
are comparab le to cluster s as tracers of thin disk char- 
acteristics, as lAfsar et al.l conclude. Taken altogether, 
these figures as well as Figure [7] show the increasing dom- 
inance of the s-process to the chemical evolution of the 
galactic thin disk. 

5. SUMMARY AND CONCLUSIONS 

The context of this work has been the recently- 
discovered trends of open cluster n eutro n c apture ele- 
ment abundances with age and R gc . D09 and IM11I have 
identified increasing [s/Fe] with decreasing age for open 
clusters that affects both heavy and light s-process ele- 
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Figure 12. [X/Fe] versus [ Eu/Fe] for Ba (top) , La (middle) and Zr ( 
horizontal branch stars from [Afsar et al.| (20123) (black crosses) in the 

ments; IY12I found such a trend only for Ba. We have also 
remarke d on the analysis of solar neighborhood dwarf 
stars by Ida Silva et al.l (|2012l ). who also reported trends 
of increasing [s/Fe] ratios with age, especially for stars 
younger than 4.5 Gyr, though we noted the contradic- 
tory _lhidingj3j3foH^ star stud- 
ies (jReddv et al.|[20?53HBensbv et al.|[2005l ). 

Here, we have presented an analysis of neutron capture 
element abundances in a sample of 19 open clusters that 
is completely independent of the open cluster studies de- 
scribed above. In this context, we summarize the results 
of this work as the following: 

• We have found a statistically significant trend 
of increasing cluster mean [Ba/Fel with decr eas- 
ing cluster ag e , in a greement with ID09I . IY12I and 
Ida Silva et all (pOll) . 

• In contrast, [La/Fe] and [Zr/Fe] abundances for our 
sample do not show trends with age. The dist ribu- 
tion of [L a/Fe] and [Zr/Fe] abundances of the [MTTl 
and IY12I samples are not inconsistent with our re- 
sults. 

• We report a marginally significant trend of in- 
creasing cluster mean [Eu/Fe] with increasing clus- 
ter age for our sample, while trends of [Ba/Eu] 
and [La/Eu] versus age are statistically signifi- 
cant. This indicates the increase in dominance of 



m). Symbols are the same as in Figure [8] along with field red 
le panels. Dotted lines represent the 1-to-l relation. 

s-process enrichment to the thin disk over time. 

• Investigations of [s/Fc] versus R gc for our clus- 
ter sample found no statistically significant trends, 
bearing in mind the limited R gc range of our sam- 
ple (^9-13 kpc). That said, clusters older than 4 
Gyr do show a marginally statistically significant 
increase in [Eu/Fe] versus R gc that the younger 
clusters do not. This may be an indication of the 
flattening of disk abundance gradient s with time 
(jFriel et al.ll200l IJacobson et al.|[201lT) . 

• The ratio of heavy to light s-process element abun- 
dances ([hs/ls]) versus [Fe/H] for our clu ster s am- 
ple is flat, consiste nt with the results of IM11I and 
IReddv et all (|2003D for thin disk dwarfs. One in- 
terpretation of this result is that the efficiency of 
the s-process in AGB stars is independent of their 
m etallicity in this me tallicity range, as suggested 
bv lReddv et all (p003l) . 

• The younger clusters in our sample show enhanced 
[s/r] ratios compared to the older open clusters, as 
one might expect based on the trends of individ- 
ual elements with age discussed above. There is 
no trend of [s/r] ratio with [Fe/H], again indicating 
that metallicity does not seem to play a crucial role 
in evolution of neutron capture elements in open 
clusters (or more explicitly, in open clusters span- 
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Figure 13. [Ba/Eu] (top panels) a nd [La/Eu] (bottom panels) versus [Fe/H] for our sample. Symbols are the same as in Figure [8] and 
again we include the RHB stars from Afsar ct al. (2012) as small crosses. 



ning the metallicity range of our sample: [Fe/H]~ 
-0.30 to +0.10). The clusters' [s/r] ratio is also 
seen to be consistent with other disk populations 
(RHB stars). 

How does one reso lve the inc onsis tenci es between the 
results of this study, ID09I/IM11I and IY121 as well as the 
inconsistencies in the behaviors of the heavy and light 
s-process elements? Are these all due to errors or a 
too- large dispersion in, e.g., our Zr and La abundances 
that mask trends that are easily seen for Ba? Given the 
strength of the Ba lines and the sensitivity of the EW 
measurements to continuum placement and/or blend- 
ing, as well as their sensitivity to uncertainties in mi- 
croturbulent velocity, we consider our La measurements 
to be more reliable. It is possible that the lack of clus- 
ters younger than 700 Myr in our sample and in IY12I 
may limit the detectability of trends with age if only 
the youngest objects have large enhancements in [s/Fe]. 
It is difficult to disentangle the impacts of different age 
ranges and systematic abundance differences on the con- 
clusions found by different studies. Clearly, it imper- 
ative to expand the analysis to include other neutron 
capture species (both s- and r-process) in not only this 
sample, but in a larger open cluster sample. We look for- 
ward to the Gaia-ESO Survey of stars in many dozens of 
open clusters that will facilitate the measure of many n- 
capture species in clu sters spanning a wide range in age 
(jGilmore et al.l [20T2T) . We hope that the existence and 



strength of any abundance trends with age so far seen 
(or not) can be verified in a larger, homogeneous study 
of open cluster chemical abundances. 
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